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A US  I'  HACT 

Progress  on  the  Distributed  Sensor  Networks  program  is  reported. 
Topics  covered  in  the  sensor  area  are  the  performance  of  acoustic 
sensors,  sizing  of  acoustic  front-end  computation  requirements, 
feasibility  of  adapting  packet  radios  to  perform  as  radar  sensors, 
and  a preliminary  investigation  of  IR  options.  In  all  cases,  the 
targets  considered  were  low-flying  aircraft  or  cruise  missiles. 
Progress  is  reported  in  the  area  of  decision  theoretic  surveillance 
space  search  methods  of  multisite  detection  and  on  specific  algo- 
rithms to  locate  targets  using  two  acoustic  sites.  Results  of  a 
preliminary  parametric  sensor  and  system  cost  analysis  are  given. 
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I.  INTRODUCTION  AND  SUMMARY 

This  is  the  second  Semiannual  Technical  Summary  (SATS)  report  of  the  Lincoln  Laboratory 
Distributed  Sensor  Networks  (DSN)  program  which  is  aimed  at  identifying  and  demonstrating  in- 
novative applications  of  new  developments  in  computer  network  theory  and  computer  science  to 
systems  employing  multiple  sensors  for  target  surveillance  and  tracking.  Such  systems  are 
made  up  of  sensors,  data  bases,  and  processors  distributed  throughout  an  area,  interconnected 
by  an  appropriate  communication  system,  and  serving  users  who  are  also  distributed  and  served 
by  the  same  communication  system.  A major  hypothesis  to  be  demonstrated  is  that  through  suit- 
able netting  and  distributed  processing,  the  information  from  individual  sensors  of  limited  ca- 
pacity and/or  range  can  be  combined  to  yield  an  effective  survivable  surveillance  system  at  fea- 
sible cost.  The  program  is  oriented  around  the  eventual  development  of  a testbed  system,  and 
demonstrations  in  the  context  of  such  a testbed.  The  testbed  problem  selected  is  the  detection 
and  tracking  of  low-flying  aircraft  and  cruise  missiles.  This  is  an  important  application  prob- 
lem since  terrain  masking  and  clutter  pose  major  difficulties  for  classical  systems  based  upon 
larger  ground-based  surveillance  radars. 

I'rogress  and  results  obtained  during  this  reporting  period  are  summarized  below  and  pre- 
sented in  greater  detail  in  subsequent  sections  of  this  report. 

A.  SENSORS 

Our  accomplishments  in  four  sensor-related  areas  may  be  summarized  as  follows. 

(1)  We  completed  further  evaluation  of  the  performance  of  acoustic  sensors, 
particularly  acoustic  arrays.  Acoustic  azimuth  measurements  from 
digital  processing  were  demonstrated  over  a complete  10-km  aircraft 
track.  This  was  done  in  the  presence  of  coherent  directional  inter- 
fering noise.  The  ability  to  discriminate  on  the  basis  of  direction  and 
frequency  content  was  thus  also  demonstrated.  We  also  verified  that 
transmission-path  variations  make  it  very  difficult  to  directly  obtain 
aircraft  elevation-angle  information  except  when  the  elevation  angle 

is  very  large.  Measurements  of  coherence  between  microphones 
displaced  by  0.1,  km  indicated  that  long  baseline  coherent  combining 
of  raw  acoustic  data  will  probably  not  be  useful. 

(2)  We  formulated  and  sized  a number  of  strawman  algorithms  for  detec- 
tion at  a small  acoustic  array.  In  general,  it  appears  that  the  front- 
end  signal  processing  and  detection  load  will  be  too  large  for  a typical 
minicomputer,  but  that  it  is  easily  within  the  capability  of  a minicom- 
puter coupled  with  a modest  amount  of  special-purpose  signal -processing 
hardware. 

(1)  We  completed  a study  to  determine  the  feasibility  of  adapting  packet 

radio  to  provide  radar-sensor  information  to  the  DSN.  A large  number 
of  options  were  considered.  In  general,  we  concluded  that  packet  radio 


1 


I 

i 


could  serve  as  a useful  sensor  for  low-flying  aircraft  or  cruise  missiles 
only  after  undergoing  very  major  changes  with  respect  to  power  levels, 
signal  waveforms,  and  antenna  directionality.  Even  then,  success  might 
lie  marginal  due  to  the  frequency  band  of  operation.  The  recommended 
approach  is  to  develop  effective  small  radars  whose  outputs  would  be 
distributed  as  required  by  colocated  packet  radios  and  interpreted  by  the 
DSN  as  a whole. 

(4)  We  investigated  the  possible  use  of  IK  sensors,  and  established  from  the 
literature  and  from  analysis  that  targets  of  interest  will  produce  detect- 
able signals  at  ranges  of  several  kilometers.  However,  a substantial 
sensor-development  effort  would  be  required  to  determine  if  an  appro- 
priate DSN  IK  sensor  could  be  developed.  Critical  engineering  factors 
are  cost,  complexity,  and  the  ability  to  operate  unattended  for  long 
periods  of  time. 

B.  MULTISITE  DETECTION 

Work  continues  on  the  general  problem  of  multisite  detection,  and  particularly  upon  the  use 
of  acoustic  azimuth  measurements. 

Experimental  software  has  been  developed  and  tested  to  investigate  the  characteristics  and 
performance  of  the  decision  theoretic  surveillance  space  search  approach  to  multisite  detection 
and  location.  A number  of  space  and  time  quantization  issues  have  been  identified  for  future 
study.  Using  this  software,  we  can  now  examine  the  capability  of  an  acoustic  system  to  function 
with  multiple  targets,  false  alarms,  and  errors  in  measurements. 

A companion  effort  to  develop  specific  algorithms  to  locate  targets  given  only  the  time 
history  of  azimuth  observations  from  two  acoustic  sites  has  been  successful.  Regions  of  poor 
location  capability  have  been  identified,  as  has  the  effect  of  aircraft  altitude.  These  results 
provide  insight  into  the  performance  of  an  acoustic  system.  The  method  may  be  used  as  an 
alternative  or  supplement  to  surveillance  space  search  - most  likely  it  is  a supplement.  This 
algorithm  is  not  intended  to  provide  optimum  locations  or  tracking.  It  makes  no  assumptions 
about  the  dynamics  or  trajectories  of  targets. 

C.  GENERAL  SYSTEM  ANALYSIS 

We  have  completed  an  analytical  study  of  sensor  cost  and  performance  trade-offs.  Factors 
considered  were  cost,  sensor  options  for  improving  signal-to-noise  ratios  (SNRs),  signal- 
propagation  losses,  and  sensor  ranges.  Within  the  context  of  the  study  and  the  models  used, 
the  optimum  choice  of  sensor  range  was  determined.  In  general,  this  study  was  intended  to 
identify  and  clarify  some  of  the  issues,  it  considered  only  a small  number  of  the  relevant  factors. 
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Work  has  continued  in  the  area  of  sensor  options  for  the  low-flying-aircraft  problem. 
Acoustic  and  rarlar  options  continue  to  be  investigated.  We  also  initiated  an  investigation  of  the 
feasibility  of  infrared  sensors. 

The  DSN  acoustic  sensor  option  remains  the  most  developed  option  at  this  time.  Additional 
data  from  our  fort  lluachuca,  Arizona  field  experiment  were  digitized  and  analyzed  to  further 
demonstrate  the  basic  capabilities  of  these  sensors  to  detect,  locate,  and  separate  multiple  tar- 
gets. Strawman  algorithms  for  coherent  acoustic  array  data  reduction  have  been  sized  to  verify 
that  the  required  computation  is  not  excessive  and  to  indicate  the  scale  of  hardware  require- 
ments for  an  acoustic  node  of  a DSN. 

Our  work  in  the  radar-sensor  area  has  intentionally  been  closely  coupled  to  the  character- 
istics of  packet  radios.  We  initiated  and  have  essentially  completed  a study  to  determine  if 
packet  radios  could  be  modified  or  evolved  into  units  which  could  furnish  sensor  as  well  as  com- 
munication functions  for  a DSN.  Such  a study  was  clearly  required  before  investigation  was  be- 
gun into  other  radar  options  that  do  not  utilize  packet  radio.  We  generally  conclude  that  it  is 
not  possible  to  evolve  packet  radio  into  a sensor  without  very  major  changes  to  the  packet  radio. 
However,  our  study  clearly  indicated  that  the-  design  of  a low-power  small-scale  DSN  radar  for 
low-flying  aircraft  can  probably  be  developed.  A separate,  unconstrained  study  of  DSN  radar 
options  is  recommended. 

The  investigation  of  DSN  fit  sensors  for  low-flying  aircraft  was  limited  to  establishing  that 
targets  do  produce  detectable  signals  at  ranges  of  several  kilometers.  An  in-depth  sensor  de- 
sign effort  would  be  required  to  determine  if  an  lit  sensor  appropriate  for  DSN  use  can  be  devel- 
oped. The  critical  issues  are  cost,  complexity,  and  ability  to  operate  unattended  for  long  peri- 
ods of  time. 

Results  of  these  sensor  investigations  are  summarized  in  the  following  sections. 

A.  ACOUSTIC  -SKNSOK  I'KRKO  KMANCK 

Initial  results  of  the  fort  lluachuca  acoustic  experiment  were  presented  in  our  first  SATS.1 
Visual  analysis  of  strip-chart  recordings  established  that  acoustic  sensors  have  sufficient  de- 
tection ranges  to  be  of  considerable  use  in  a DSN  for  low-flying-aircraft  surveillance.  Process- 
ing a limited  amount  of  digital  data  indicated  that  useful  acoustic  azimuth  measurements  could 
be  obtained  for  acoustic  source-receiver  distances  out  to  at  least  several  kilometers.  Since 
that  time,  additional  digital  data  covering  complete  30-km  test  tracks  have  been  obtained.  New 
processing  results  using  the  new  data  as  well  as  some  of  the  previously  available  data  are  pre- 
sented in  thiR  section. 

In  general,  we  have  concerned  ourselves  with  the  performance  of  an  acoustic-only  or  an 
acoustic -aided  DSN  for  the  surveillance  of  multiple  low-flying  targets  in  the  presence  of  noise 
due  to  either  natural  or  artificial  sources.  To  a large  extent,  acoustic-sensor  performance 
depends  upon  the  ability  of  array  processing  to  separate  incoming  acoustic  waves  according  to 
their  arrival  azimuth  and  phase  velocity  and  upon  spectral  analysis  to  differentiate  and  charac- 
terize the  nature  of  the  source  producing  the  acoustic  energy. 


TABLE  II- 1 


ARRAY  PROCESSING  RESULTS  FOR  THE  COMPLETE  A -7  TRACK 

Acoustic 

Acoustic 

Phase 

Acoustic 

Radar 

Range 

Velocity 

Azimuth 

Azimuth 

(km ) 

(m/sec ) 

(deg) 

(deg) 

-16.0 

321 

270 

264 

-11.3 

371 

249 

264 

-6.8 

318 

270 

264 

-4.0 

329 

261 

262 

-0.7 

429 

239 

240 

0.7 

371 

1)1 

120 

1.4 

364 

104 

107 

2.0 

331 

96 

102 

2.6 

329 

99 

99 

3.2 

323 

96 

98 

3.9 

359 

96 

96 

4.5 

359 

96 

95 

5.2 

341 

97 

95 

5.8 

341 

97 

94 

6.4 

361 

90 

93 

7.  1 

36) 

90 

91 

7.5 

342 

90 

91 

8.7 

342 

90 

91 

9.9 

350 

90 

91 

11.1 

328 

90 

91 

12.3 

417 

90 

91 

13.2 

375 

90 

90 

14.4 

i 

382 

80 

90 

I.  Small-Aperture-Array  Performance 

Table  II - 1 shows  digital-data-analysis  results  for  a complete  jet-aircraft  track.  The  basic 
analysis  tool  was  frequency-wavenumber  processing.  The  microphone  array  consisted  of  9 sen- 
sors in  a 5-m  aperture.  The  jet  aircraft  flew  approximately  200  m off  the  ground  over  the  en- 
tire path.  The  data  were  processed  to  obtain  azimuth  measurements  at  approximately  5-sec 
intervals.  Each  processed  data  segment  consisted  of  1 sec  of  signal  which  was  analyzed  with 
1-llz  resolution  at  10-llz  intervals  over  the  band  80  to  130  Hz.  From  the  table,  it  can  be  seen 
that  the  aircraft  was  successfully  tracked  from  16  km  inbound  to  14.4  km  outbound.  These 
ranges  are  not  necessarily  the  limiting  tracking  ranges,  but  correspond  to  the  earliest  and  lat- 
est data  recorded  in  the  field.  Note  that,  even  though  the  data  were  analyzed  in  approximately 
equal  time  intervals,  the  range  bins  are  not  equispaced.  This  results  from  the  moving  source 
effect  discussed  in  our  last  SATS.*  The  acoustic  range  was  derived  from  the  recorded  radar 
tracking  data. 

figure  Il-A-1  plots  the  observed  acoustic  azimuth  against  the  acoustic  range  calculated 
from  the  radar  tracking  information.  The  calculated  acoustic  azimuth  of  the  aircraft  is  indi- 
cated by  the  solid  line.  While  a minor  amount  of  scattering  is  present,  the  overall  performance 
of  the  processing  is  good.  In  fact,  the  acoustic  azimuths  fall  within  the  precision  of  the  radar 
data  over  the  complete  track.  The  largest  azimuthal  deviations  occur  on  the  extremes  of  the 
path.  In  these  regions  the  signal  is  more  than  40  dB  down  from  its  peak  amplitudes  near  the 
closest  point  of  approach  (CPA).  In  this  situation,  we  are  probably  suffering  from  the  limita- 
tions of  the  analog  recording  method  used  to  obtain  the  digital  data.  Nevertheless,  azimuth 
measurements  could  be  made,  even  though  (as  shown  below)  a strong  interfering  noise  back- 
ground, which  is  several  decibels  above  the  aircraft  noise,  is  present. 

Figure  II-A-2  shows  the  signals  recorded  across  the  array  at  the  time  the  acoustic  wave 
from  the  jet  at  a 16-km  distance  inbound  would  arrive.  Clearly,  a large  coherent  signal  is  prop- 
agating across  the  array  and  good  array  processing  results  should  be  expected.  Figure  Il-A-3 
shows  the  spectra  of  these  signals  in  the  80-  to  130-Hz  band.  (Unfortunately,  the  system  noise 
outside  this  band  was  significant  for  the  signal  levels  of  interest.)  The  narrowband  peaks  in  the 
spectra  indicate  that  the  signal  is  dominated  by  a piston-type  source  rather  than  the  broadband 
signal  one  would  associate  with  a jet  source.  Each  of  the  peak  frequencies  was  investigated  in 
the  wavenumber  domain  and  found  to  originate  from  a constant  direction.  The  results  at  100  Hz 
are  shown  in  Fig.  II-A-4,  where  the  arrows  indicate  the  direction  to  the  A-7  and  the  direction 
to  the  generator  used  in  the  experiment.  The  sound  in  this  interval  is  dominated  by  the  gener- 
ator noise.  Since  we  expect  the  jet  to  have  acoustic  energy  distributed  over  the  entire  range, 
the  wavenumber  spectrum  at  a point  between  the  harmonic-generator  noise  peaks  was  calculated. 
Figure  II-A-5  shows  the  result  at  90  Hz,  where  the  signal  power  is  down  by  about  10  dB  from  the 
generator  peaks.  As  expected,  the  A-7  is  detected  with  some  interference  from  the  generator 
sound  field.  In  the  region  up  to— 6 km  and  beyond  14  km,  the  dominant  signal  was  generator 
noise.  Fortunately,  because  of  the  spectral  nature  of  the  noise  and  the  signal,  successful  pro- 
cessing at  a SNR  of  at  least  —10  dB  was  obtained.  In  general,  because  source  type  has  a great 
influence  on  signal  spectra,  these  results  should  be  regarded  as  typical  and  not  unusual.  Im- 
proved array  design  and  higher  quality  data  collection  will  result  in  even  better  processing  gains. 

In  theory,  the  elevation  angle  of  an  aircraft  should  be  indicated  by  the  phase  velocity  of  the 
acoustic  wave  crossing  the  array.  The  observed  phase -velocity  values  for  the  above  experiment 


are  given  in  Table  11-1.  Despite  the  fact  that  the  array  resolution  in  phase  velocity  is  about  a 
meter/second,  the  observations  fluctuate  several  tens  of  meters/second  with  respect  to  their 
predicted  values  based  upon  radar  tracking  and  known  aircraft  altitude.  This  scatter  corre- 
sponds to  a standard  deviation  in  elevation  angle  of  more  than  10°  over  the  entire  track.  Such 
a scatter,  probably  caused  by  the  velocity  microstructure  of  the  atmosphere,  is  too  large  for 
precision  elevation-tracking  purposes.  Useful  elevation-angle  information  can  only  be  obtained 
when  the  true  elevation  angle  is  large  relative  to  the  scatter. 

2.  Uarger-Aperture-Array  Options 

One  question  we  have  not  previously  addressed  is  whether  coherent  combining  of  raw  acous- 
tic data  (as  opposed  to  combining  single-site  azimuth  or  frequency  measurements)  from  sites 
separated  by  several  kilometers  is  possible.  To  some  extent  this  is  still  an  open  question,  but 
we  do  have  some  results  which  indicate  that  it  may  not  be  possible  or  practical.  Fortunately, 
the  ability  to  do  such  long-range  coherent  combining  is  not  crucial  for  the  functioning  of  a DSN. 
f urther  research  in  this  area  might  be  done  in  parallel  with  other  future  experimental  and  test- 
bed efforts. 

Using  jet-aircraft  (A-7)  and  helicopter  (UH-1)  data,  we  looked  at  iong-range  coherence  by 
computing  the  cross -cor relation  between  central  array  elements  and  the  CUD  microphone  offset 
from  the  array  by  600  m (see  Ref.  1 for  a complete  description  of  the  experiment).  The  600-m 
distance  is  not  as  large  as  one  would  like  to  answer  the  questions,  but  it  was  the  largest  which 
could  be  obtained  during  the  experiments.  Figure  II-A-6  shows  the  cross-correlation  of  the 
GRA  and  GRC  elements  of  the  central  array.  These  microphones  are  located  only  0.5  m apart. 
The  signal  was  from  a UH-1  helicopter.  As  expected,  the  peak  correlation  is  high  (0.9)  and  the 
envelope  of  the  cross-correlation  function  decays  very  slowly  due  to  the  narrowband  nature  of  the 
sound.  The  cross-correlation  between  the  GRA  and  GRD  instruments  is  shown  in  Fig.  II-A-7; 
these  instruments  are  located  600  m apart.  Time  differences  of  arrival  corresponding  to  wave- 
front  group  velocities  equal  to  the  speed  of  sound,  c,  and  equal  to  600  m/sec  are  indicated  in 
the  figure.  While,  as  one  would  predict,  the  envelope  of  the  function  peaks  in  this  time-lag  re- 
gion, the  envelope  drops  too  slowly  to  be  certain  which  peak  to  select.  The  dominant  energy  in 
the  signal  is  about  40  Hz  so  that  the  wavelength  is  about  10  m.  The  maxi  mum  correlation  is  0.4. 
Thus,  over  approximately  60  wavelengths  the  signal  peak  coherence  drops  by  more  than  a fac- 
tor of  2.  In  a poorer  signal-to-noise  environment,  the  resolution  would  be  even  worse.  We  are 
led  to  the  result  then  that,  over  several  tens  of  wavelengths  for  narrowband  signals,  acoustic 
path  characteristics  degrade  signals  to  the  point  where  coherent  signal  processing  would  pro- 
vide minimal  signal  information.  Figures  II-A-8  and  -9  show  the  same  computation  for  the  A-7. 
Jet  acoustic  sources  are  broadband,  and  the  improvement  in  resolution  is  obvious  even  though 
the  overall  peak  coherence  is  smaller  (0.3  as  opposed  to  0.4  for  the  UH-1)  than  for  the  narrow- 
band  source.  These  preliminary  results  indicate  that  coherence  of  both  broadband  and  narrow- 
band  sources  will  drop  at  sensor  separations  of  a kilometer  or  more  to  the  point  that  coherent 
processing  of  raw  data  will  not  be  possible  at  all.  On  a preliminary  basis,  we  conclude  that  the 
systematic  combination  of  multiple  small  arrays  will  not  be  improved  with  long-range  coherence 
of  the  actual  acoustic  signals.  Fortunately,  it  should  be  possible  to  design  and  demonstrate  an 
acoustic  DSN  without  requiring  such  long  baseline  coherent  processing.  Further  investigations 
of  long  baseline  coherent  processing  could  be  a parallel  longer-term  research  activity.  Coher- 
ent processing  of  spectra  (as  opposed  to  time  series)  from  widely  separated  sites  is  another  area 
which  could  be  studied  in  the  future. 
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Fig.  II-A-1.  Azimuth  detections  over  the  complete  A-7  track.  200  m altitude. 
Negative  ranges  are  inbound. 
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Fig.  II-A-2.  Acoustic  signal  recorded 
when  A-7  was  16  km  inbound.  Note 
high  signal  coherence  level  across  ar- 
ray. This  narrowband  coherent  signal 
is  noise  from  a nearby  generator.  The 
A-7  signal  is  not  visible. 


Fig.  II-A-3.  Spectra  on  a linear  scale 
of  signals  shown  in  Fig.  II-A-2.  Nar- 
rowband peaks  indicate  signals  are 
dominated  by  piston-type  source  rather 
than  jet  source.  P.  aks  are  about  10  dB 
above  background. 
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n.  ACOUSTIC  COMUUTATIONAI.  SI/.ING 

Thi-  current  concept  for  an  acoustic  sensor  in  a DSN  is  that  each  sensor  site  will  include  a 
small  array  of  microphones.  Many  such  sensor  sites  will  be  distributed  over  the  area  of  inter- 
est. Kach  small  array  will  consist  of  several  microphones  situated  within  an  aperture  of  a few 
meters.  A considerable  amount  of  processing  must  be  performed  on  the  raw  acoustic  data  at 
each  small  array  to  detect  and  calculate  an  initial  estimate  of  the  apparent  azimuth  of  incoming 
acoustic  signals.  It  is  important  to  estimate  the  computational  requirements  (speed  and  mem- 
ory) of  this  front-end  processing  to  see  if  a hardware  configuration  that  fulfills  the  requirements 
is  feasible,  and,  if  it  is,  to  estimate  the  size,  nature,  and  complexity  of  the  hardware  that  might 
be  deployed  in  the  field  with  the  acoustic  arrays.  Multisite  detection  and  tracking,  signature 
analyses,  or  knowledge-based  interpretation  of  data  are  not  considered  in  the  sizing  results 
repo  rted. 


1.  Strawman  Algorithms 

Three  illustrative  strawman  algorithms  were  chosen  for  the  sizing.  The  computation  re- 
quirements of  these  algorithms  are  typical  of  those  that  would  be  used  in  a front-erid  processor 
for  acoustic  sensors.  The  three  algorithms  have  been  named  (a)  Incoherent  Frequency-Domain 
Detection,  (b)  Coherent  Frequency-Domain  Detection  or  Frequency-Domain  Beamforming.  and 
(c)  Time- Domain  Beamforming.  The  first  algorithm  only  detects  the  presence  of  targets,  while 
the  latter  two  algorithms  detect  and  also  calculate  a rough  azimuth  for  each  target.  Frequency- 
Domain  Beamforming  is  a variation  of  the  frequency-wavenumber  processing  we  have  used  to 
analyze  small  array  data  from  our  Fort  Huachuca  experiment.  The  strawman  algorithms  are 
briefly  defined  below.  In  each  case  the  input  data  are  N-point  sampled  waveforms  from  each  of 
M sensors: 

{x^n)}  for  i = t,2,...,M  and  n=  1,2 N 

a.  Incoherent  Frequency-Domain  Detection  (Target  Detection  Only) 

(1)  Fourier  transform  each  of  the  waveforms: 

Xjfn)  =>  X.(n) 

(2)  Calculate  the  magnitude  squared  of  all  the  positive-frequency  components  of  the 
spectrum: 

| Xjfn)  | 2 for  i = 1,2,...  ,M  and  n = 1,  2 N/2  . 

(3)  Average  the  spectrum  components  for  each  positive  frequency  over  the  M sensors 
of  the  array: 

M 

Y(n)  = ^ Yj  l^ln)!2  • 

i=  1 

(4)  Smooth  the  Y values  using  an  s-point  rectangular  kernel  and  resample  in  frequency 
at  a lower  rate  to  obtain  Z(q),  q = 1,2,...  ,N/2r  where  r is  the  decimation  factor 
for  resampling  the  smoothed  spectrum. 
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(5)  Update  the  long-term  time  average  of  power  at  each  resampled  frequency  using  a 
single-pole  recursive  filter  which  results  in  exponential  weighting  of  past  values 
in  the  average; 

W(q)  I W(q)  l '/.( q)  for  q 1,2 N/2 r 

(6)  Calculate  a detection  statistic: 

N/2r 

I)  /.(q)/W(q)  . 

q=l 

When  the  detection  statistic  exceeds  a threshold,  a target  is  declared  to  be  present. 
This  particular  detection  statistic  is  just  the  sum  of  detection  statistics  for  each  of 
the  frequencies.  The  development  of  actual  practical  algorithms  to  declare  detec- 
tions based  upon  the  individual  detection  statistics  is  an  area  for  serious  future 
research.  We  expect  that  actual  detection  algorithms  may  lie  come  quite  complex 
and  involve  many  heuristics.  This  also  is  true  of  the  other  detection  algorithms 
described  in  the  following  sections. 
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b.  Coherent  Frequency-Domain  Detection  (Frequency- Domain 
Beamforming) 

(1)  Fourier  transform  the  waveforms: 
x.(n)  =S-  X.(n) 


(2)  Calculate  a power  estimate  1*  for  each  azimuth  k and  for  each  frequency  j under 

JK 

consideration.  This  is  done  using  precalculated  steering  vectors  C^.  Steering 
vectors  are  complex  vectors  whose  elements  represent  the  phase  delay  re- 
quired to  align  the  signal  for  a particular  frequency  j,  azimuth  k,  and  sensor  i 
with  the  signal  from  other  sensors.  The  calculation  is: 


M 

'V  ^ Gjk,xiW  • 

1=1 

The  number  of  azimuths  is  a function  of  frequency  and  is  selected  to  give  overlap 

in  array  beams  at  the  3-dB  points.  The  number  of  azimuths  used  for  calculation 

was  2jrfR/C  where  f is  frequency,  K is  the  diameter  of  the  acoustic  array,  and 

C is  the  speed  of  sound.  Only  the  Q frequencies  in  a preselected  band  from  f . 

to  f are  considered.  The  number  of  azimuths  of  interest  at  the  highest  fre- 
max 

quency  is  V. 

(3)  Smooth  power  estimates  in  frequency  using  an  s-point  rectangular  smoothing  ker- 
nel. Kesample  the  smoothed  power  estimates  saving  every  r**1  frequency  point. 

This  smoothing  is  done  with  fixed  azimuth.  Since  the  number  of  azimuths  is  a func- 
tion of  frequency,  we  cannot  exactly  do  the  frequency  smoothing  for  fixed  azimuth 
but  we  can  do  it  approximately  provided  the  smoothing  kernel  is  not  too  broad  in 
frequency.  The  resulting  quantities  are  where  q goes  from  1 to  Q/r  and  for 

each  q,  m goes  from  1 to  the  number  of  azimuths  considered  at  the  frequency 
corresponding  to  q. 
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(4)  Take  long-term  averages  of  the  smoothed  and  resampled  power  estimates  using  a 
single-pole  recursive  filter: 

W =cW  \ 
qm  qm  qm 

(5)  Calculate  detection  statistics  corresponding  to  each  of  the  V azimuths  of  interest 
at  the  highest  frequency  of  interest: 

q=Q/r 

v i 

q-1 

where  m(q,  v)  indicates  the  azimuth  index  which  at  frequency  q (more  precisely 
the  frequency  corresponding  to  index  q)  corresponds  to  the  azimuth  closest  to  azi- 
muth v (more  precisely  the  azimuth  corresponding  to  index  v).  When  a detection 
statistic  exceeds  a threshold,  a target  is  declared  at  the  v1*1  azimuth.  As  noted 
previously,  the  definition  and  use  of  actual  detection  statistics  is  a matter  for 
serious  future  research. 

c.  Time-Uomain  Beamforming 

(1)  Beamform  in  the  time  domain;  that  is,  add  the  waveforms  from  the  M sensors 

with  appropriate  delays  to  form  V (see  Table  II— 2)  beams  of  N points  each.  Each 
beam  represents  the  acoustic  signals  propagating  through  the  array  from  a differ- 
ent azimuth: 

Xj(n)  , i = 1.  Z M , n = 1,  Z N 

=t-  yv(n)  , v = 1,  Z V , n = 1 , Z N . 


^qm(q,v) 

qm(q.v) 


TABLE  11-2 

SENSOR  CONFIGURATION  AND  PROCESSING  REQUIREMENTS 
FOR  STRAWMAN- ALGORITHM  SIZING 

M = 7 sensors 

N = 512  waveform  samples  from  each  sensor  per  second 

R = 3 m (diameter  of  sensor  orray) 

f = 200  Hz  (maximum  frequency  of  interest) 
max 

f . = 5 Hz  (minimum  frequency  of  interest) 
min  ' 

Q = 196  frequencies  of  interest  in  algorithm  b 

r = 10  (frequency  resampling  interval) 

s=10  (length  of  frequency  smoothing  kernel) 

V = 20  beams  (the  number  of  resolvable  beams  at  f which  is  2*f  R/c  = 11.4, 
where  c is  the  speed  of  sound)  mOX  mOX 
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(2)  Fourier  transform  each  beam: 


yv(n)  =i>  Yv(n)  . 

(3)  Calculate  the  magnitude  squared  of  the  positive  frequency  components  of  the  spec- 
trum for  each  beam  as  was  done  in  step  (2)  of  algorithm  a.  for  each  of  the  individ- 
ual channels. 

(4)  Smooth  and  resample  each  beam  as  done  for  the  array  average  powers  in  step  (4) 
of  algorithm  a. 

(5)  Update  the  long-term  averages  for  each  beam  using  the  same  single-pole  recur- 
sive filter  as  in  step  (5)  of  algorithm  a.  for  each  beam  spectrum. 

(6)  Calculate  detection  statistics  for  each  beam  which  are  equivalent  in  computation 

to  that  used  in  step  (6)  of  algorithm  a.  When  a detection  statistic  exceeds  a thresh- 
old, a target  is  declared  to  be  present  at  the  azimuth  of  that  beam. 

To  size  and  compare  the  computational  requirements  of  the  three  algorithms,  the  sensor 
configuration  and  processing  requirements  described  in  Table  II-2  are  utilized.  The  assump- 
tions made  about  the  processor  are; 

• All  computations  are  done  with  16-bit  integer  arithmetic. 

• An  addition  or  subtraction  requires  one  operation  or  instruction  epoch. 

• A multiplication  requires  four  of  these  epochs. 

• A division  requires  twelve  epochs. 

• Input  data  are  double-buffered.  Thus,  the  time  required  for  reading  the 
data  is  assumed  to  overlap  with  processing  time. 

• All  calculations  are  done  on  the  acoustic  data  for  all  the  sensors  in  the 
array  once  every  second. 

2.  Results 

The  major  results  (the  number  of  operations  per  second  and  storage  requirements  for  pro- 
gram and  data)  as  well  as  the  breakdown  of  the  various  computational  tasks  are  shown  in 
Table  11-3.  The  formula  for  the  number  of  operations  to  compute  the  fast  Fourier  transform 
assumes  that  an  N/2-point  complex  FFT  can  be  done  on  N real  data  points.  It  is  also  assumed 
that  a complex  butterfly  requires  24  operation  epochs  (4  multiplies  and  8 additions)  and  that 
15-percent  overhead  is  required  for  bit-reversal,  shuffling  of  the  data,  and  other  details. 

The  first  algorithm,  target  detection  only,  can  be  performed  within  the  given  time  con- 
straints with  a medium-sized  minicomputer  such  as  the  PDP  11/34.  However,  the  second  two 
algorithms  would  seriously  tax  or  exceed  such  a minicomputer's  computational  capacity.  Thus, 
a medium-sized  minicomputer  alone  would  not  be  sufficient  for  a front-end  processor  in  an 
acoustic  array. 

A significant  fact  is  that  in  all  three  algorithms,  the  major  fraction  of  computational  time 
is  spent  performing  the  FFT:  92  percent  of  the  instructions  in  algorithm  a.,  54  percent  in  b., 
and  80  percent  in  c.  A relatively  fast  array  processor  (with  approximate  maximum  FFT  com- 
putational capacity  of  a 256  complex-point  FFT  in  5 msec)  could  significantly  alleviate  the  pro- 
cessing load  on  a front-end  central  computer. 
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Am  example  of  a system  that  could  easily  handle  the  processing  requirements  of  all  the  algo- 
rithms is  shown  in  Pig.  tt -11  - 1 . This  system  does  not  necessarily  represent  the  configuration 
that  would  he  used  at  a node  in  a DSN  in  the  field,  hut  rather  represents  a research  vehicle  that 
could  he  used  to  further  develop  and  refine  acoustic-sensor  detection  and  location  algorithms. 

The  implementation  includes  a packet  radio  for  communication,  a minicomputer  (the  1’DD  11/34), 
disk  storage  (for  program  and  data  storage),  a signal  processor,  and  analog-to-digital  convert- 
ers. The  configuration  has  significant  reserve  capacity  such  that  a wide  variety  oi  computational 
options  could  lie  explored. 

The  processor  shown  in  l ig.  II-H-1  is  the  Lincoln  Digital  Signal  Processor  (I. DSP)  which 
is  capable  of  20-million  instructions  per  second,  hut  has  a data-storage  capacity  oi  only 
4,090  10-hit  words  and  a program -storage  capacity  of  2,040  10-hit  words.  However,  using 
current  memory  technologies,  an  I, DSP-type  processor  can  he  designed  with  more  data  storage 
(for  example.  05,530  10-bit  words)  with  a small  sacrifice  in  speed  - a factor  of  2 or  3.  While 
acoustic  target  detection  algorithms  alone  do  not  require  such  a sophisticated  fast  processor, 
one  might  lie  appropriate  to  perform  processing  for  several  different  types  of  sensors  in  one 
area  (for  example,  small  radars,  infrared,  and  acoustic  sensors).  There  are  other  less- 
sophisticated  and  less-expensive  processors,  including  those  using  recently  developed  charge- 
coupled-device  (CCD)  technologies,  that  are  currently  available  and  can  be  investigated  for  this 
application. 


C.  1'ACKF.T  RADIO  AS  A RADAR  SENSOR 


In  order  for  a Racket  Radio  Net  (PRN)  to  provide  digital  communication  among  a set  of 
widely  scattered  ground  terminals  (some  pairs  of  which  are  beyond  line-of-sight  (LOS)],  the 
packet  radios  would  have  to  be  deployed  fairly  uniformly  over  the  area,  with  spacings  limited 
to  keep  nearest  neighbors  within  LOS.  Depending  on  terrain,  foliage,  and  siting,  radio  spacings 
would  vary  from  a few  to  perhaps  as  much  as  20  km.  A relatively  uniform  grid  of  radios  would 
provide  multiple  connectivity  for  each  radio  and  allow  for  a reliable  adaptable  communication 
system.  A similar  geographic  deployment  of  radar  sensors  would  also  provide  good  low-altitude 
coverage  for  detection  of  low-flying  aircraft  and  cruise  missiles.  We  report  here  the  results  of 
a study  to  assess  the  feasibility  of  modifying  or  evolving  the  packet  radios  to  function  as  radar 
sensors,  as  well  as  communication  terminals. 

1.  Distributed  Radar  Sensor  Net  Options 

The  packet  radio  communication  waveforms  themselves  are  well  suited  to  accurate  timc- 
of-arrival  (TOA)  measurement,  which  is  needed  for  bit  synchronization  to  decode  a packet. 

This  TOA  measurement  capability  will  also  be  used  in  a later  version  of  a PRN  to  provide  for 
automatic  position  location  of  terminals.  The  TOA  measurement  made  on  a packet  signal  re- 
flected from  an  airborne  target  would  provide  a measure  of  target  range.  The  range  resolution 
provided  by  the  Upgraded  Packet  Radio  (UPR)  waveform  would  be  on  the  order  of  3 m.  The  sim- 
plest radar  sensor  might  thus  be  an  omnidirectional  I PR  which  measures  only  range  delay  of 
targets.  Monostatic  (transmission  and  reception  of  target  signals  at  the  same  site)  or  bistatic 
radar  modes  would  both  be  useful.  Range-delay  measurements  from  three  omni-sensors  would 
be  required  to  obtain  a target- location  estimate  in  three  dimensions.  This  implies  a need  for 
overlapping  coverage  areas  of  adjacent  sensors.  Of  course,  a directional  antenna  could  also 
be  employed  for  the  radar  function  to  provide  range  and  azimuth  information  on  targets  as  a 
conventional  radar  does,  but  this  option  entails  a more-significant  evolution  of  the  UPR.  Both 
monostatic  and  bistatic  ranging  sensors  have  been  considered  in  this  study,  as  well  as  omnidi- 
rectional and  narrow-beam  antennas. 

Although  the  UPR  signal  is  modulated  by  data  and  a bandspreading  key,  it  is  possible,  at 
least  in  principle,  to  process  a data  packet  with  arbitrary  modulation  and  detect  the  presence 
of  a doppler-shifted  signal  due  to  a moving  target.  One  way  to  do  this  would  be  to  square  the 
received  signals,  producing  phase -coherent  pulses  at  twice  the  "mark'1  and  "space"  frequencies 
of  the  original  MSK  signal.  Filtering  would  have  to  be  employed  to  separate  the  direct  path  sig- 
nal from  the  reflected  (doppler-shifted)  signal  for  target  detection.  This  mode  of  operation  is 
more  analogous  to  that  of  a CW  radar  measuring  doppler  information  only.  A packet  transmis- 
sion of  10  msec  would  allow  a maximum  (theoretical)  doppler  resolution  of  100  Hz,  which  at  the 
UPR  frequency  of  1.8  GHz  corresponds  to  a doppler  velocity  of  16.7  m/sec.  Practically  achiev- 
able doppler  resolution  would  be  poorer.  The  best  doppler  resolution  must  also  be  a small  frac- 
tion of  the  maximum  doppler  caused  by  a moving  target,  implying  that  useful  doppler  informa- 
tion may  be  achievable  with  targets  moving  200  m/sec  or  more  (390  knots). 

As  an  alternative  to  processing  the  packet-data -transmission  format,  a special  CW  signal 
could  also  be  transmitted  somewhere  in  the  UPR  band  to  be  used  as  a CW  radar  signal.  The 
simplest  form  of  doppler-only  sensor  is  again  the  omnidirectional  sensor  which  could  measure 
a doppler-vs-time  profile  for  a target.  For  targets  with  constant  but  unknown  velocity,  three 
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such  doppler-vs-ttme  profiles  are  needed  [or  a target-track  estimate  in  three  dimensions.  As 
in  the  consideration  of  ranging  sensors,  both  monostatic  and  histatic  dopplor-only  sensors  are 
considered  as  well  as  omnidirectional  and  narrow-beam  antennas. 

In  summary,  the  set  of  options  considered  in  this  study  may  be  described  as  follows: 

Measurements  ltangc  vs  doppler-only 

Mode  Monostatic  vs  bistatic 

Antenna  Omnidirectional  vs  narrow  beam 

Waveform  Ill'll  waveform  vs  special  waveforms 

Considering  the  bistatic  mode  antenna  options  to  give  four  separate  cases,  a total  of  24  sensor 
options  have  been  considered  in  this  study. 

2,  Target  Model 

The  targets  of  interest  to  a distributed  U I'd  sensor  net  are  low-flying  aircraft  and  cruise 
missiles,  l or  baseline  calculations  the  fighter-aircraft  target  is  considered  first,  since  cruise 
missiles  are  smaller  and  will  be  harder  to  detect,  f igure  Il-C-1  shows  the  monostatic  radar 
cross  section  (HCS)  of  a jet  fighter  aircraft  in  a straight  and  level  attitude  with  the  radar  in  the 
azimuth  plane.  The  nose-on  HCS  is  generally  the  lowest,  presenting  the  most  difficult  detection 
situation.  HCS  measurements  indicate  that  jet  fighter  aircraft  typically  present  between  5-  and 
10-dBsm  KCS  in  the  nose-on  sector  at  the  Ill'll  frequency.  The  target  model  adopted  for  base- 
line calculations  in  this  study  is  10  dBsm,  which  is  considerably  larger  than  would  be  rxp>  ted 
for  nose-on  views  of  cruise  missiles. 

fighter  aircraft  attempting  to  penetrate  enemy  air  defenses  at  low  altitudes  would  probably 
not  exceed  speeds  of  Mach  1 or  accelerations  of  4 to  6 Cl's. 

3.  Baseline  SNK  Calculations 

The  following  power  budget  assumes  Ill'll  parameters  and  a 10-dBsm  target  10  km  from 
transmitter  and  receiver  (l.H-Gllz  frequency). 

Transmitter  power  110  W)  = 10  dBW 

Transmit  antenna  gain  = 9 dll 

Path  loss  outbound  (10  km)  = —117.5  dB 
Target  gain  (10  dBsm)  = 30.5  dB 

Path  loss  inbound  (10  km)  = —117.5  dll 
deceive  antenna  gain  = 9 dBl 

deceived  signal  power  = —170.5  dBW 

for  an  8-dB  noise-figure  receiver,  the  noise  power  density  is  Nq  = —190  dBW/Hz.  Thus,  for 
a single  10-psec  UPR  pulse,  the  SNd  is 

f /N  =- 170.5  dBW  - 50  dB  sec  + 190  dBW/llz  = -24.5  dB 
p o 

Since  an  SNd  of  at  least  +10  dB  is  required  for  detection,  the  above  calculation  indicates  that  a 
3 5-dB  improvement  in  SNd  is  required  beyond  the  omnidirection  UPR  parameters  to  detect  a 
10-dBsm  target  at  10  km  with  a single  10-psec  pulse. 
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Monostatic  and  bistatic  (overage  areas  consistent  with  the  10-kxn  target  range  from  trans- 
mitter and  receiver  are  shown  in  Fig.  II-C-2.  Overlapping  coverage  areas  require  sensor  spac- 
ings  of  10  km  or  less. 

The  SNR  improvement  options  are; 

(a)  Increased  transmitter  power, 

(b)  Lower  noise  receiver  front  end, 

(c)  l.onger  coherent  integration  time, 

(d)  Antenna  gain,  transmit  and/or  receive. 

The  transmitter  power  obtainable  from  a solid-state  transmitter  at  1.8  GHz  is  currently  limited 
to  a few  hundred  watts  (average).  The  receiver  noise  figure  could  be  reduced  to  about  3 dB  with 
a new  !•' FT  front  end.  Thus,  an  improved  solid-state  transmitter  (say  300  W)  and  front  end 
could  provide  20  dH  of  the  required  3 5-dB  SNit  improvement.  In  an  omni-sensor,  the  remain- 
ing 13  dH  of  improvement  would  have  to  be  obtained  by  extending  the  coherent  integration  time 
to  at  least  300  nsec. 

A modest  directional-beam  antenna  (10  azimuth  beamwidth)  could  provide  a mainbeam  gain 
of  about  10  dH  more  than  the  omnidirectional  UPR  antenna.  Directional  antennas  on  both  trans- 
mit and  receive  could  thus  provide  20  dH  of  SNR  improvement  which,  together  with  the  improved 
transmitter  and  receiver,  would  allow  detection  of  the  desired  target  with  a single  10-psec  pulse. 

Thus,  there  are  means  available  to  provide  adequate  SNR  for  the  aircraft  targets  of  inter- 
est at  reasonable  detection  ranges.  However,  there  are  other  complicating  factors  which  must 
be  considered  beyond  the  SNR,  some  of  which  are: 

(a)  Coherent  integration  of  doppler-shifted  signals, 

(b)  Ground  clutter, 

(c)  Jamming, 

Id)  Compatibility  of  directional  antennas  with  UPR  channel  access  protocol. 

4.  Jamming 

Because  the  received  power  from  a target  varies  as  the  fourth  power  of  the  target  range, 
a radar  is  much  more  vulnerable  to  jamming  than  a communication  link.  Figure  II-C-3  shows 
a plot  of  jamming-to-signal  ratio  (JSR)  for  an  omnidirectional  sensor  receiver  vs  target  range 
for  various  ratios  of  jammer  to  radar  transmitter  effective  radiated  power  (ERP).  This  plot 
shows  that  even  for  modest  jammer  ERPs  (relative  to  UPR),  the  JSR  requires  more  improve- 
ment than  the  SNR  for  the  UPR  used  as  a radar  sensor.  It  is  important  to  note  that  a single 
airborne  jammer  at  a reasonable  standoff  range  would  be  able  to  simultaneously  disrupt  many 
omnidirectional  radar  sensors  spread  out  over  a wide  area.  We  concluded  that  all  sensor  op- 
tions which  include  omnidirectional  receivers  are  thus  infeasible  because  of  the  vulnerability 
to  a single  airborne  jammer. 

3.  Transmitter/Receiver  Isolation 

A radar  which  transmits  a signal  of  sufficiently  long  duration  to  overlap  the  received  echo 
signal  from  a target  requires  isolation  of  the  transmitter  from  the  receiver.  A ranging  sensor 
using  a long-duration  coded  waveform  (like  the  UPR  packet  waveform)  would  require  suppres- 
sion of  the  transmitted  signal  sufficiently  to  reduce  the  correlation  sidelobes  of  the  transmitted 
signal  below  that  of  the  correlation  peak  of  the  target  echo.  A CW  doppler  sensor  would  require 
transmitter  suppression  sufficient  to  reduce  the  spectral  sidelobes  of  the  transmitted  signal 
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bo  low  that  of  the  spec  tral  peak  of  the  target  echo,  l or  a monostatic  radar  using  a single  antenna 
for  transmit  and  receive,  the  transmitter-power-to-reccived-echo-power  ratio  (T HR)  at  the  an- 
tenna port  would  be  two  dll  for  a 10-dllsm  target  at  10  km.  In  this  case,  it  is  not  feasible  to  ob- 
tain enough  transmitter  suppression  to  detect  the  target  signal.  1 his  rules  out  the  monostatic 
radar  options  using  nonpulsed  signals  (i.e.,  those  options  using  'he  I I’lt  packet  waveiorm  or  a 
CW  signal  for  floppier  estimation). 

A significant  amount  of  transmitter-receiver  isolation  can  be  realized  by  employing  sepa- 
rate transmit  and  receive  antennas,  which  really  amounts  to  a bistatic  radar  mode.  I he  case 
of  a monostatic  radar  using  separate  hut  closely  spaced  antennas  has  therefore  not  been  treated 
as  a separate  case. 

An  omni-bistatic  radar  with  transmitter  and  receiver  separation  of  10  km  will  experience 
a THU  of  80  dll.  The  correlation  sidelobes  of  a ranging  sensor  using  a nonpulsed  coded  wave- 
form such  as  the  packet  waveform  could  only  be  a maximum  of  00  dll  below  the  correlation  peak 
for  an  ideal  coded  waveform  extending  to  the  maximum  coherent-integration  limit  oi  10  msec. 
Practical  levels  of  correlation  sidelobes  are  more  on  the  order  of  30  to  40  dll.  putting  the  cor- 
relation peak  of  a radar  echo  40  to  SO  dll  below  the  sidelobes  of  the  direct  signal.  The  only 
practical  bistatic  ranging  radar  option  that  uses  an  extended  waveform  would  be  one  employing 
directional  antennas  on  both  transmitter  and  receiver.  The  direct-path  signal  would  then  be 
suppressed  by  the  main-beam-to-sidelobo-gain  ratio  of  both  antennas,  which  could  be  on  the 
order  of  30  dH  per  antenna.  The  correlation  sidelobes  of  the  direct  signal  need  then  be  only 
30  dl3  below  the  correlation  peak  to  allow  detection  of  the  desired  target. 

A CW  doppler-only  bistatic  radar  will  liave  to  utilize  the  maximum  possible  coherent  inte- 
gration time  (10  msec)  in  order  to  achieve  the  desired  doppler  resolution  (iOO  Hz).  1 his  implies 
that  the  direct-path  signal  must  be  suppressed  sufficiently  to  reduce  the  near-in  spectral  side- 
lobes of  the  direct  signal  below  the  spectral  peak  of  a target  echo.  For  the  minimum  separation 
in  frequency  of  100  liz  between  the  direct-path  signal  and  a doppler-shifted  echo,  it  would  be  dif- 
ficult but  possible  to  suppress  the  direct  signal  relative  to  the  echo  as  much  as  60  dB  under  ideal 
conditions.  Thus,  the  only  bistatic  doppler-only  radar  option  that  appears  practically  feasible 
is  also  one  employing  directional  antennas  on  both  transmitter  and  receiver. 

In  summary,  transmitter-receiver  isolation  for  radar  modes  using  long-duration  signals 
which  would  overlap  received  echoes  is  a very  significant  factor  which  eliminates  many  of  the 
radar  options. 

6.  Summary  and  Conclusions 

Table  11-4  lists  the  twelve  ranging-sensor  options  that  were  considered  in  this  study.  The 
table  shows  that  jamming  vulnerability  and  the  transmitter-receiver  isolation  requirements  taken 
together  eliminate  all  but  four  of  the  options.  Table  £1-5  lists  the  twelve  CW  doppler-only  sensor 
options  studied.  Jamming  eliminates  all  options  using  the  UPR  waveform  because  the  operation 
of  all  radars  on  a common  carrier  frequency  increases  the  jamming  vulnerability  of  the  set  of 
sensors.  The  transmitter-receiver  isolation  requirements  also  eliminate  an  additional  number 
of  options,  resulting  in  only  one  surviving  option.  The  surviving  options  are  listed  in  Table  11-6 
in  order  of  increasing  difficulty  of  implementation.  It  appears  that  the  changes  required  to 
evolve  packet  radios  into  practical  radar  sensors  would  be  very  major.  Moreover,  the  detection 
of  low-flying  aircraft  and  cruise  missiles  is  a difficult  sensor  problem,  even  with  no  constraints 
upon  the  sensor  design.  Rather  than  try  to  evolve  packet  radios  into  a radar  sensor,  we  recom- 
mend that  a small  radar  sensor  be  developed  to  be  colocated  with  packet  radios  in  a DSN. 


TABLE  11-4 


RANGING-SENSOR  OPTIONS  TREF 
(Feasible  Options  Marked  With  • ) 

Mode 

Transmit 

Antenna 

Receive 

Antenna 

Waveform 

ECM 

Vu  Inerabi  lity 

Direct-Path 
i Suppression 
(dB) 

1 

Required  SNR 
Improvement 
(dB) 

UPR 

Severe 

190 

35 

i 

Special 

Severe 

Pulsed 

35 

Monostatic 

1 

Narrow  Beam  Narrow  Beam 

UPR 

Operable 

: 

190 

15 

Special 

Operable 

Pulsed 

’5  . 

t ' ' 

UPR 

Severe 

90 

35 

Omni 

Special 

Severe 

Pulsed 

35 

Bistatic 

Narrow  Beam 
1 

UPR 

Moderate 

60 

25 

Special 

Moderate 

Pulsed 

25  . 

UPR 

Severe 

60 

25 

Narrow  Beam 

umm 



Special 

Severe 

Pulsed 

25 

Narrow  Beam 

UPR 

Operable 

30 

'5  . 

Special 

Operable 

Pulsed 

15  . 





DOPPLER-ONLY 

(Feasible 

TABLE  11-5 

(CW)  SENSOR  OPTIONS  TREE 

Option  Marked  With  • ) 

Mode 

Transrr.it 

Antenna 

Receive 

Antenna 

Waveform 

ECM 

Vulnerabi  lity 

— 

Direct-Path 

Suppression 

(dB) 

Required  SNR 
Improvement 
(dB) 

Monostotic 


Narrow  Beam  Narrow  Bean 


Operable 


Narrow  Beam 


Bistatic 


Moderate 


Narrow  Beam 


Narrow  Beam 


Operab  le 


D.  SKAKCH-MODi:  HASSIVK  IK  FOK  I ,OW- 1 T.Y  INC  AlltCKAFT 

In  this  section,  we  examine  the  use  of  lit  detectors  in  the  context  of  a DSN  whose  objective 
is  to  detect  and  track  low-flying  aircraft.  In  the  spirit  of  our  previous  acoustic  investigations, 
we  restrict  the  sensor  to  be  passive  in  its  basic  mode  of  operation  and  to  be  capable  of  unattended 
operation.  The  problem  naturally  divides  itself  into  three  areas:  (1)  sources  and  background 
noise  that  might  obscure  or  lie  misidentified  as  sources,  (2)  source-receiver  path  effects,  and 
H)  receiver  characteristics.  The  basic  questions  we  wish  to  answei^are  what  can  be  detected 
and  at  what  distance  can  detection  be  made.  Here,  we  will  try  to  utYswer  these  questions  by  ex- 
amining the  basic  engineering  principles  of  Ilt  in  the  context  of  a specific  example. 


equation 


Detection  range,  target,  path,  and  detector  characteristics’  are  related  through  the  range 
2 


WA  t = 4111/2 

nil2  a HDD*  t NT 
o 


= Source  power  iri  the  band  of  interest, 

It  Detection  range, 
t&  = 1,oh8  in  the  atmosphere, 

SI  Total  solid  angle  scanned, 

D = Diameter  of  the  optics, 

D*  = Detectivity  of  the  sensor, 

t - 1,088  due  to  optics, 
o 

N = Number  of  detectors, 

T = Time  to  scan  If,  and 

SNIt  = Signal-to-noise  ratio  required  for  a given  probability 
of  detection. 

The  term  on  the  left-hand  side  is  the  power  flux  that  arrives  at  the  optics  due  to  a source  at 
distance  It.  The  left-hand  side  describes  the  power  flux  required  at  the  aperture  to  match  the 
detector-optics  syBtem  noise  multiplied  by  the  SNIt  that  will  produce  the  desired  confidence 
level  for  detection. 

Airborne  targets  emit  energy  ubove  normal  background  levels  from  the  heated  skin,  hot 
engine  parts,  and  hot  exhaust  gases.  Generally,  each  of  these  sources  has  different  radiation 
characteristics,  any  one  of  which  might  be  best  for  detection  purposes  depending  on  the  aspect 
and  dynamics  of  the  target,  the  construction  of  the  vehicle,  or  atmospheric  conditions.  The 
total  power  flux  produced  by  an  object, 

W = <oT4 

is  proportional  to  the  fourth  power  of  its  absolute  temperature  T and  its  emissivity  t,  where  a 
(the  Stefan-Holtzmann  constant)  is  S.f>7  * 10"12  W/(cm2-K4).  At  supersonic  speeds  the  hot  skin 
may  dominate  the  radiation,  whereas  at  lower  speeds  the  relative  size  of  the  exhaust  plume  and 
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hot  engine  parts,  their  temperatures,  and  emissivities  determine  the  proportion  of  radiation 
from  each.  The  distribution  of  radiant  power  at  various  wavelengths  A is  given  by  Planck's 
equation 

WA  = t^CjdA/^A1’  {exp[Cz/(AT)|  - 1}) 
where  is  the  ennssivity  at  the  wavelength  A,  and 

C'j  = 3.75  < 1012  W-cm2 

C'2  = 1.438  104  pm-K 

l or  constant  emissivity.  between  room  temperature  (300  K)  and  the  melting  point  of  steel 
(2900  K),  the  peak  of  the  distribution  decreases  from  about  10  to  1 pm  spanning  the  infrared 
spectrum.  Since  the  atmosphere  severely  attenuates  certain  energy  bands  and  since  the  com- 
position of  the  exhaust  as  well  as  its  temperature  and  shape  determine  its  radiation,  an  exact 
statement  of  source  radiation  is  difficult  to  make.  Generally  speaking,  the  radiation  pattern  is 
the  superposition  of  the  three  effects  combined  with  the  shielding  effect  of  the  aircraft  itself. 

The  result  is  commonly  an  order-of-magnitude  difference  in  radiated  power  at  different  azimuths 
and  inclination  angles.  As  an  example,  jet  aircraft  generally  radiate  between  100  and  1000  W 
total  power  with  plume  temperatures  of  about  1000  K (see  Kef.  3).  If  a source  is  emitting  100  W 
total,  then  in  a 1 -pm  band  about  2 pm,  the  power  emitted  is  about  10  W as  can  be  derived  from 
the  above  equations  assuming  a 1000  K source.  At  a distance  of  10  km  then,  the  power  flux  is 

t W./irll2  = t * 3 * 10 "41  W/cm2 

cl  A U 

4 

The  2-pm  band  was  chosen  as  an  example  since  plume  radiation  is  strong  in  the  band  and  for 
the  detector  and  path  considerations  that  are  explained  below. 

In  the  lower  atmosphere,  water  and  CO.,  severely  attenuate  the  propagation  of  IK  energy. 
However,  conditions  exist  that  allow  certain  spectral  bands  of  energy  to  be  transmitted  over 
large  distances.  In  the  IK  region,  the  bands  of  transmittance  are  8 to  14  pm,  3 to  5 pm,  and 
0.7  to  2.5  pm  (see  Kef.  5).  In  each  of  these  bands,  attenuation  still  exists  to  the  degree  of  about 
0.5  dli/km  on  the  best  of  days  up  to  many  dll/km  during  rain,  fog,  or  haze  conditions.  Specif- 
ically, near  2 pm.  when  the  visibility  is  about  10  km  the  attenuation  is  about  4 dll/km.  At  25-km 
visibility,  the  attenuation  is  about  2 dll/km.  Following  upon  the  example  above,  4-dB/km  atten- 
uation at  10  km  leads  to  a power  flux  of 

W.  la  ? 

t — = 3 * 10‘”  W/cm  . 
a *K2 

It  remains  then  to  examine  detectors  suitable  for  DSN  applications. 

The  IK  receiver  system  consists  of  the  optics  for  collecting  and  focusing  incoming  target 
photons  onto  a detector  that  counts  them.  The  instantaneous  field  of  view  (1FOV)  of  the  system 
may  be  estimated  by  considering  the  background  noise  characteristics  that  can  lead  to  misidenti- 
fying  false  targets  (such  as  cloud-top  reflections)  with  real  targets.  For  instance,  in  the  2-pm 
band  the  worst  background  level  from  sun  reflection  at  the  earth's  surface  will  be  about 
1.5  x 104  W/cm2-sr  (see  Kef.  2).  Since  reflections  are  diffuse  compared  with  real  targets. 
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discrimination  between  them  will  be  possible  il'  the  IKOV  is  small  compared  with  the  ratio  of  the 
target  power  flux  at  the  aperture  to  the  reflection  flux  per  steradian.  f or  the  example,  we  find 
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3 ■ 10 -1 3 W/cni2  , .<) 

IFOV  = i j ^ 10  sr 

1.5  10  W/cm  -sr 

or  about  45  prad  on  a side,  which  is  about  the  resolving  power  of  5-(  m circular  optics/'  Thus, 
a telescope  with  5-cm  optics  or  larger  is  required.  The  focal  depth  of  the  telescope  will  depend 
on  detector  size  requirements.  If  a search  is  conducted  at  a single  elevation  at  all  azimuths, 
then  the  total  solid  angle  scanned  will  be 

Si  = 2jt  45  10'^'  sr  = 2.8  10~4  sr 


Of  course,  many  elevations  could  be  scanned  utilizing  a vertical  array  of  sensors.  Say,  for 
tracking  purposes,  that  we  wish  to  scan  Si  every  second  (an  aircraft  traveling  near  the  speed 
of  sound  moves  about  300  m in  1 sec).  Thus,  the  time  t to  dwell  at  each  point  in  the  field  of 
view  would  be 


t-  -2  ‘ 10  >r.  l sec  * 7 


2.8  10  sr 


psec 


IK  detectors  are  characterized  by  their  detectivity  D*  which  is  extremely  dependent  on  de- 
tector temperature,  time  constant,  and  spectral  response.  Given  the  parameters  above,  the 
range  equation  can  be  rearranged  to  study  detector  characteristics  trade-offs.  Assuming  an 

optics  loss  of  0.5  for  a 10-cm  aperture  and  a SNR  of  5 (this  assures  99.9-percent  probability  of 
2 

detection  in  a single  scan  ),  we  find 


D*  (cm/W-sec1/2)  N3' 2 = 3.7 


1033  cm/(W-sec3/2) 


Table  II -7  shows  two  possible  detectors  that  might  be  used  in  a system.  An  array  of  four 
cooled  indium  antimonide  detectors  would  suffice,  while  it  would  take  60  x S0D  lead  selenide  de- 
tectors to  accomplish  the  task.  l.ead  sulfide  is  a common  room-temperature  detector  with  a 
large  detectivity,  but  it  unfortunately  has  a 2 50 -psec  rise  time  and  is  thus  unsuitable  since  the 
dwell  time  required  for  our  example  is  only  7 psec. 


The  above  analysis  is  not  an  attempt  at  a sensor  design.  What  is  clear  from  the  above  is 
that  low-flying  aircraft  produce  a sufficient  photon  flux  at  many  kilometer  distances  to  be  de- 
tected, and  therefore  lit  sensors  should  be  further  considered.  However,  the  above  also  sug- 
gests that  the  system  performance  requirements  for  IK  sensors  may  not  be  met  by  simple  or 
inexpensive  sensors  due  to  cooling  or  very  large  detector  array  fabrication  requirements. 
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III.  Ml  I.TISiTE  DETECTION 


In  the  area  of  multisite  detection,  the  utilization  of  acoustic  azimuth  measurements  is 
emphasized.  A DSN  composed  of  multiple  small  acoustic  arrays  for  detection  and  tracking  of 
low-flying  aircraft  seems  most  appropriate  for  the  first  stages  of  a testhed  system.  It  is  thus 
clear  that  more  understanding  of  the  multisite  utilization  of  such  acoustic  data  is  called  for. 

Two  complementary  efforts  were  undertaken,  and  progress  is  summarized  in  the  following 
sections.  The  first  was  to  develop  the  surveillance  space  search  approach  suggested  in  our 
previous  SATS.  In  that  effort,  some  basic  experimental  software  has  been  developed  and  some 
initial  evaluation  has  been  done  for  the  case  of  acoustic  sensors.  Several  issues  needing  fur- 
ther study  have  been  identified.  Computational  sizing  has  been  deferred  until  a planned  major 
restructuring  of  the  algorithms  is  done.  The  second  effort  was  to  develop  and  evaluate  algo- 
rithms to  directly  locate  aircraft  using  azimuth  observations  from  only  two  acoustic  sites 
measuring  azimuth.  For  this  it  was  assumed  that  the  sensor  site  could  associate  a sequence 
of  azimuth  observations  with  a single  target  so  that  it  actually  delivered  a time  series  of  azimuth 
observations  already  associated  with  each  other. 

A.  DECISION  THEORETIC  SURVEILLANCE  SPACE  SEARCH 

MUt/riSITE  DETECTION 

In  this  section,  we  report  progress  on  general  algorithms  that  combine  data  from  several 
sensor  sites  to  detect  and  locate  targets  even  when  data  from  a single  sensor  may  not  be  suffi- 
cient to  locate  a target.  This  work  has  been  done  in  the  context  of  acoustic  sensors  which  can 
measure  the  acoustic  azimuth  of  targets.  For  such  acoustic  sensors  the  azimuth  from  which 
sound  is  arriving  may  not  correspond  to  the  azimuth  at  which  the  target  is  located.  This  is  a 
factor  which  complicates  the  problem  and  must  be  addressed.  The  approach  we  evolved  can  be 
adapted  to  a wide  variety  of  sensors  without  propagation  delay  problems,  and  can  also  handle 
multiple-sensor  types  without  added  difficulty. 

The  surveillance  space  search  approach  to  the  problem  is  to  divide  the  physical  space  of 
interest  into  three-dimensional  cells  (the  space  covered  by  the  cells  may  or  may  not  include 
the  location  of  sensors),  and  then  search  the  space  cell -by-cell  for  the  targets.  For  the  acous- 
tic problem  of  immediate  concern,  we  assume  that  the  sensors  periodically  make  azimuth 
measurements,  and  record  the  time  of  the  measurement  along  with  the  azimuth.  We  also  as- 
sume that  the  sensors  are  synchronized  (to  tenths  of  seconds)  with  each  other  and  the  multisite- 
detection  processor. 

The  search  for  targets  is  made  through  each  cell  in  the  target  space  at  periodic  intervals 
called  surveillance  times.  The  search  technique  is  as  follows.  For  each  cell,  the  hypothesis 
is  made  that  there  is  a sound-emitting  target  within  the  boundaries  of  the  cell  at  the  surveillance 
time.  Given  a cell,  the  data  from  each  sensor  covering  the  surveillance  space  are  considered 
in  turn.  The  azimuth  that  each  sensor  would  measure  if  there  were  indeed  such  a target  in  the 
cell  at  the  surveillance  time,  and  the  time  of  arrival  of  this  sound  at  the  sensor  are  calculated. 
The  azimuth  and  time  are  actually  ranges  of  azimuth  and  time  because  the  cell  has  finite  dimen- 
sions. This  range  of  azimuth  is  then  compared  with  those  data  that  the  sensor  actually  measured 
during  the  predicted  time  range.  If  there  is  an  agreement,  a count  is  listed  in  this  cell  for  the 
particular  sensor  at  the  surveillance  time.  Any  additional  data  that  the  sensor  may  have,  such 
as  if  this  measurement  is  consistent  with  others  comprising  the  track  of  one  particular  target. 
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i ;m  also  be  ini  billed  with  the  cell  corroboration.  When  all  the  cells  have  been  searched,  targets 
are  declared  in  those  cells  with  corroborations  from  many  sensors.  The  number  of  corrobora- 
tions is  an  ad  hoc  decision  statistic  selected  for  initial  studies.  When  the  search  is  done  at 
several  successive  surveillance  times,  targets  can  be  tracked  as  they  move  from  cell-to-cell. 

A significant  advantage  to  this  detection  approach  is  that  it  is  not  necessary  to  make  any 
assumptions  about  the  dynamics  (the  velocity  or  trajectory)  of  the  target.  One  only  assumes  a 
position  m a cell.  A disadvantage  of  such  a multisite-detection  scheme  for  acoustic  sensors  is 
that  the  searches  must  he  done  at  a delay  to  real  time.  Since  targets  are  hypothesized  in  each 
cell  at  the  surveillance  time,  one  must  account  for  the  acoustic  travel  time  of  the  targets'  sound 
from  the  cell  to  the  sensor.  Thus,  the  surveillance  time  or  time  of  the  search  must  be  done  at 
a delay  equal  to  the  (ravel  time  from  a cell  to  the  most  distant  sensor  from  which  one  desires 
azimuthal  data. 

1,  Software 

A simulation  program  of  the  surveillance  space  search  algorithm  has  been  written  in  the 
C programming  language  to  run  on  a PI.iT  I I.  The  program  consists  of  three  major  sections 
initialization,  search,  and  report. 

In  the  initialization  routines,  the  position  of  sensors  and  the  initial  positions  and  velocities 
of  the  targets  are  used  to  calculate  the  acoustic  azimuths  that  each  sensor  measures  as  the 
targets  pass  through  the  area  of  interest.  Only  targets  with  constant  velocity  have  been  used. 

A rectangular,  box-shaped  grid  of  cubed -shaped  cells  is  also  set  up  over  the  area  and  altitudes 
of  interest.  As  a result  of  these  routines,  each  sensor  has  a table  that  includes  azimuth  mea- 
surements taken  over  the  time  interval  in  which  the  targets  passed  through  the  area  of  interest 
and  the  time  for  each  measurement.  In  the  current  simulation,  it  is  assumed  that  all  sensors 
make  measurements  simultaneously,  once  a second. 

The  details  of  the  search  routines  can  be  explained  with  the  aid  of  Fig.  Ill  - A - 1 (a -b) . As 
mentioned  above,  surveillance  space  search  is  done  after  a delay  from  the  actual  real  time. 

The  sensor  data  that  are  used  to  corroborate  the  hypotheses  that  targets  are  in  the  various 
cells  include  all  those  measurements  made  from  the  time  of  the  search  to  the  current  real  time, 
as  indicated  in  Fig.  Ifl-A-lfb).  During  the  search,  each  cell  in  the  space  is  checked.  Once  a 
cell  is  chosen,  data  from  each  sensor  are  then  considered  in  turn. 

For  each  cell-sensor  combination,  two  times,  t . and  t (those  times  required  for 

m i n m qx 

acoustic  signals  to  travel  from  the  closest  and  farthest  boundaries  of  the  cell  to  the  sensor), 
are  first  calculated  as  shown  in  Fig.  Ml-A-l(a).  The  two  times  and  the  geometrical  relationship 
between  the  sensor  and  the  surveillance  cell  provide  a mapping  of  the  surveillance  space  cell 
into  a data  space  cell  in  the  sensor's  measurement  space  as  shown  in  the  figure.  If  measured 
azimuths  are  found  in  the  data  space  cell,  they  are  said  to  corroborate  the  hypothesis  that  a 
target  is  in  the  surveillance  cell  at  the  surveillance  time,  and  this  agreement  is  noted  in  the 
surveillance  cell  along  with  the  sensor  number  and  the  time  of  the  measurement(s).  In  the 
particular  example  shown  in  Fig.  II I - A - 1 (a-b),  there  are  two  azimuth  measurements  in  one  of 
the  data  cells  because  the  target  remained  within  the  surveillance  cell  during  two  measurement 
periods  of  the  sensor.  Another  surveillance  space  cell  mapping  is  shown  in  the  figure.  No 
measurements  appear  in  this  second  data  space  cell.  When  checking  if  the  azimuth  measure- 
ments fall  in  the  data  space  cell,  the  current  implementation  uses  an  error  margin  of  plus  or 
minus  one-half  degree. 
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After  all  fells  arc  searched,  a list  of  corroborations  can  be  printed  out  or  the  entire  (trod 
of  cells  can  be  displayed  with  corroborations  noted,  as  is  shown  in  Figs.  Ill- A~2  through  lll-A-H 
in  the  sections  which  follow.  'I  hose  cells  with  corroborations  from  all  or  nearly  all  the  nearby 
sensors  are  considered  to  very  likely  have  one  or  more  targets  in  them. 

2.  Target  l)ete<  tion  and  I 01  .die.. 

Some  sample  scenarios  were  run  through  the  simulation.  The  area  of  interest  consists  of 
a box,  5-km  square  and  0.5  km  high.  The  box  is  quantized  into  100  cells,  each  a 0.5-km  cube. 
Five  acoustic  sensors  are  distributed  over  the  area,  all  are  at  zero  altitude.  Four  sensors  are 
at  the  corners  of  the  1-km  square  area  of  interest  and  the  fifth  is  at  the  center. 

a.  Static  Target 

To  more  easily  understand  how  the  algorithm  works,  a surveillance  space  search  was  done 
for  a static  target  at  zero  altitude.  The  search  was  done  after  a 20-sec  delay,  that  is,  after 
the  sensors  had  recorded  azimuthal  data  for  20  sec.  The  corroborations  of  predicted  azimuths 
and  sensor  data  are  shown  in  Fig.  lII-A-2(a).  A corroboration  that  a target  is  in  the  cell  is 
designated  simply  by  the  number  of  the  sensor.  As  was  discussed  above,  those  cells  with  more 
than  one  corroboration  from  the  same  sensor  simply  contain  more  than  one  cell-sensor  distance, 
with  an  acoustic  travel  time  equal  to  the  times  that  the  sensors  made  measurements.  As  shown 
more  clearly  in  Fig.  III-A-2(b),  the  cells  with  corroborative  data  from  one  sensor  form  a straight 
line  from  the  sensor  through  the  target.  In  the  figure,  straight  lines  were  drawn  from  the  sen- 
sors through  the  target.  The  corroborations  farther  from  the  sensor  represent  later  measure- 
ments times.  The  sensor  measures  a constant  azimuth  for  20  sec.  Thus,  the  sensor  has  an 
azimuth  that  agrees  with  that  predicted  by  cells  along  the  straight  line  through  the  target  and 
sensor  for  acoustic  travel  times  from  0 to  20  sec. 

The  cell  in  which  all  the  straight  lines  intersect  contains  corroborations  from  all  the  sen- 
sors and,  indeed,  also  contains  the  target.  For  this  example,  this  cell  has  corroborations  from 
ail  5 sensors  a total  of  12.  Other  cells  have  corroborations  from  a maximum  of  3 different 
sensors  - a maximum  total  of  4. 

b.  Moving  Target 

Detection  and  location  of  a target  that  is  moving  are  more  difficult.  As  the  sensors  make 
sequential  azimuth  measurements,  the  target  is  moving  and  its  acoustic  azimuth  is  changing. 

The  resulting  line  of  cells  containing  corroborations  is  not  straight,  as  is  shown  in  Fig.  III-A-3. 

In  the  figure,  the  celts  with  one  or  more  corroborations  by  sensor  0 for  possible  target  locations 
at  time  t 0,  using  all  measurements  available  for  the  interval  10  to  F-  20,  have  been  marked 
with  solid  circles.  The  actual  target  is  traveling  in  the  direction  shown  at  a speed  of  174.9  m/sec 
(M  0.53)  at  zero  altitude.  The  target  is  in  the  cell  marked  with  an  X at  t 0.  The  acoustic 
azimuths  that  sensor  0 measures  at  t - 0 and  20  sec  are  also  shown  in  the  figure.  A curve  was 
drawn  by  eye  through  the  cells  with  corroborations. 

As  with  the  static  target,  the  cell  in  which  corroborations  exist  from  all  sensors  is  hypoth- 
esized to  contain  the  target.  Figures  MI-A-4(a)  through  (f)  show  the  results  of  searches  at 
surveillance  times,  t = 0,  4,  K,  12,  10,  and  20  sec,  respectively,  as  real  time  progresses 
from  t 20  to  40  sec.  Each  surveillance  search  is  done  with  an  interval  of  20  sec  of  sensor 
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data.  To  avoid  cluttered  figures,  only  one  corroboration  per  sensor  is  indicated  in  a cell.  The 
cell  m which  the  target  was  located  at  the  search  time  (marked  with  an  X in  the  figure)  was 
corroborated  by  all  S sensors  at  these  times. 

The  scores  by  which  the  target  was  located  are  presented  in  Table  I f I - 1 . The  two  numbers 
shown  are  the  number  of  corroborations  in  the  cell  with  the  target  and  the  maximum  number  of 
such  corroborations  in  other  cells  in  the  grid. 

for  tins  particular  example,  the  cell  in  which  the  target  was  traveling  could  he  identified 
easily  by  selecting  those  cells  with  corroborations  from  all  sensors.  These  cells  could  also 
be  identified  by  selecting  those  cells  with  the  maximum  total  number  of  corroborations  {this 
technu|ue  is  a little  weak  for  t 12  [Fig.  III-A-4(d)|  where  the  target  is  very  close  to  a cell 
boundary}.  As  will  be  seen  below,  there  are  examples  in  which  neither  of  these  counting  tech- 
niques alone  can  select  target  cells,  but  both  can  be  used  to  select  those  cells  containing  targets. 

c.  Two  Moving  Targets 

Is  this  search  technique  able  to  distinguish  among  several  targets  traveling  through  the  area 
of  interest?  The  scenario  of  Figs.  IH-A-4(a)  through  (f)  was  augmented  with  a second,  faster 
target  traveling  at  279. 5 m/sec  (M  - 0.84),  and  zero  altitude.  The  search  algorithm  was 
able  to  identify  definitively  two  cells  with  targets  for  4 out  of  the  5 surveillance  space  searches. 
The  cells  in  which  the  targets  were  located  were  not  always  corroborated  by  ail  8 sensors,  as 
shown  in  Table  I II -2  and  Figs.  Iir-A-S(a~b)  through  fff-A-8(a-b),  where  the  trajectories  of  the 
two  targets  are  shown  as  lines.  The  locations  of  the  targets  at  the  search  time  are  marked  with 
an  X.  However,  by  considering  both  the  number  of  sensors  corroborating  each  cell  and  the 
total  number  of  corroborations  for  each  cell,  both  target-containing  cells  were  located  in  the 
four  searches.  At  t - 8 sec  fFig.  Ill- A-6(a-b|  I,  using  the  above-mentioned  criteria,  the  faster 
target  is  not  correctly  located  and  a target  is  located  erroneously  in  the  cells  that  are  shaded 
in  the  figure.  Using  the  results  of  the  searches  at  t 0,  4,  12,  and  16  sec,  the  assumptions 
that  the  targets  are  moving  at  constant  velocity,  and  some  rudimentary  tracking,  one  can  cor- 
rectly locate  the  faster  target  at  t 8 sec.  Perhaps  surveillance  searches  done  more  often 
could  also  resolve  this  problem. 

Again,  for  this  example  the  surveillance  space  search  can  detect  and  locate  two  cells  with 
targets  in  the  area  of  interest.  The  detection  and  location  appear  more  difficult  with  two  tar- 
gets, and  the  question  arises  as  to  how  often  should  the  surveillance  space  search  he  conducted. 
The  nature  of  the  tracking  process  to  he  done  at  a higher  level  is  also  brought  into  consideration. 

3.  Discussion 

The  results  described  above  comprise  only  an  initial  study  of  the  surveillance  space  search 
multisite  detection  and  location  algorithm.  There  are  a large  number  of  issues  that  have  been 
raised  in  implementing  the  algorithm. 

As  presently  implemented,  the  search  algorithm  is  not  efficient.  All  the  cells  are  searched, 
regardless  of  the  chance  of  finding  targets  in  them.  The  above  examples  use  a rather  coarse 
grid  - 0.5-km  cubes  - over  a small  area.  If  a larger  area  and  altitude  were  more  finely  quan- 
tized, it  could  easily  require  the  task  of  searching  more  than  one  million  cells  for  one  surveil- 
lance space  search.  To  keep  track  of  several  targets,  the  search  might  have  to  be  repeated 
more  often.  Alternately,  using  a finer  grid  (smaller  cell  dimensions)  might  reduce  the  necessity 
of  searching  often.  Clearly,  the  computational  load  increases  rapidly  with  the  precision  with 
which  one  wishes  to  follow  the  targets  through  the  area  of  interest. 
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Other  quantization  effects  are  apparent.  Given  that  acoustic  sensors  make  measurements 
once  a second,  a faster  target  will  result  in  fewer  measurements  as  it  travels  through  the 
surveillance  space.  Fewer  measurements  result  in  fewer  corroborations  with  the  azimuths 
predicted  by  the  cell -sensor  geometries.  One  possible  solution  is  to  interpolate  the  target's 
trajectory  - that  is.  interpolate  values  between  measurements  and  use  the  interpolated  values 
in  the  search  algorithm.  However,  the  interpolation  procedures  add  to  the  computational  load, 
also,  if  there  are  several  targets  (several  azimuth  measurements  for  one  given  time),  it  is  not 
always  clear  which  azimuths  should  be  used  for  interpolation.  An  alternative  solution  is  to  have 
the  acoustic  sensors  generate  azimuth  measurements  more  often;  however,  this  solution  results 
in  a greater  computational  load  for  the  front-end  processors. 

Further  study  is  required  to  better  understand  the  properties  of  the  search  algorithm  and 
the  utility  of  its  output.  For  example,  the  above  results  have  not  included  the  effects  of  false 
alarms  or  noisy  azimuth  measurements.  Would  several  false  azimuths  cause  the  detection 
and  location  algorithm  to  fail?  Could  a sophisticated  tracking  algorithm  use  the  results  from 
several  inconclusive  surveillance  space  searches  and  track  a target9 

Another  issue  is  target  altitude.  Targets  will  be  at  different  and  perhaps  varying  altitudes. 
Can  results  with  satisfactory  accuracy  be  obtained  from  a search  at  only  one  altitude,  or  must 
the  search  be  carried  out  through  all  the  cells  of  several  layers  of  surveillance  space,  increasing 
the  computational  load?  Work  reported  in  Sec.  B on  p.-iO  has  direct  bearing  upon  these  questions. 


DATA  USED  IN  SEARCH  AT  TIME  t„ 


Fig.  HI- A-t . Multisite  detection:  surveillance  space  search, 
(a)  Surveillance  space,  (b)  azimuthal  data  measured  at  sensor. 
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(b)  Surveillance  time  = 4 sec;  real 
time  - 24  sec. 


(c)  Surveillance  time  = 8 sec;  real 
time  = 28  sec. 


Fig.  III-A-4(a-f).  Number  of  corroborating  sensors.  Target  moving 
through  surveillance  space  at  174.9  m/sec.  Target  location  at  sur- 
veillance time  marked  by  an  X. 
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Fig.  III-A-S.  Two  moving  targets:  surveillance  time  4 sec  (real 
time  24  sec),  (a)  Number  of  corroborating  sensors;  (b)  total  num- 
ber of  corroborations. 
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l*'ig.  III-A-6.  Two  moving  targets  surveillance  time  H sec  (real 
time  28  sec),  (a)  Number  of  corroborating  sensors,  (b)  total  num- 
ber of  corroborations. 


Fig.  lll-A-8.  Two  moving  targets  surveillance  time  t6  sec  (real 
time  sec),  (a)  Number  of  corroborating  sensors,  (b)  total  num- 
ber of  corroborations. 


14.  TWO-SKNSOK  TAKGF.T  LOCATION 


A second  effort  to  develop  multisite  acoustic  location  algorithms  has  been  under  way  cori- 
current  with  the  investigation  of  decision  theoretic  surveillance  space  search  alternatives.  Our 
immediate  goal  was  to  show  that  an  aircraft  can  he  located  iri  the  horizontal  plane  using  data 
from  only  two  acoustic  sensors.  F.ach  sensor  is  a microphone  array  whose  outputs  are  of  the 
form  "sound  is  coming  from  azimuth  (horizontal  direction)  d>  at  time  t."  We  felt  that  under- 
standing how  to  do  this  might  furnish  insight  into  the  capabilities  of  acoustic  sensors  and  that 
the  algorithms  developed  might  actually  be  of  use  in  a !>SN.  The  two-sensor  method  described 
below  can  he  used  as  an  alternative  to,  supplement  to,  or  theoretical  background  for  surveillance 
space  searching. 

1.  Locating  an  Aircraft 

Suppose  we  have  two  sensors,  Sj  and  S^.  located  10  km  apart  and  both  at  the  same  height, 
which  is  zero  in  our  coordinate  system.  Suppose  a plane  is  flying  around  somewhere,  and  the 
azimuth-vs-time  curves  for  the  two  sensors  are  as  shown  in  Fig.  Ill-B-1. 

First,  we  want  to  find  the  plane  at  time  t.  assuming  (with  no  particular  reasoni  that  the 
plane  is  at  height  z 0.  Consider  just  the  S,  data.  The  sound  emitted  by  the  plane  at  time  t 
must  arrive  at  S.  at  some  time  t1  t.  Let  <t>  1 be  the  azimuth  reading  measured  for  time  t( . 
Then,  at  time  t the  plane  is  distance 

( t j - t)  • (speed  of  sound) 
from  Sj , or  at  horizontal  distance 

from  S(.  Given  t^ , we  have  found  <t>t  and  r^ , which  specify  a point  in  the  horizontal  plane. 

The  set  of  all  pairs  (rt,  <p t)  given  by  all  t(  > t generates  a curve  of  possible  positions  in  the 
horizontal  plane  for  the  airplane  at  time  t.  We  call  this  latter  curve  the  possible  position  curve 
for  Sj,  t,  and  z 0.  Figure  IJI-B-2  shows  such  curves  for  t 40,  -26,  22.  18,  14,  and 

10  sec,  z 0,  and  both  and  S^.  Note  that  the  possible  position  curve  for  a sensor  originates 

at  that  sensor,  which  is  the  point  on  the  curve  corresponding  to  t(  t,  unless  the  available  mea- 
sured azimuth  curve  begins  after  time  t. 

At  time  t,  the  plane  must  be  at  the  intersection  of  the  S,  and  S2  curves  for  given  t and  z, 
assuming  that  the  plane  is  actually  at  height  z.  Thus,  for  each  t we  get  a plane  location. 
Connecting  these  locations  with  straight  lines  gives  the  computed  plane  path  displayed  in 
Fig.  Ill- 14-2  which  is  just  about  equal  to  the  actual  plane  path,  so  the  actual  path  is  not  displayed. 
The  plane  is  actually  flying  a straight  line  at  Mach  0.9  and  height  0.1  km.  The  difference  be- 
tween the  assumed  height  z 0 and  actual  height  0.1  km  results  in  negligible  error. 

The  possible  position  curves  for  t -22  and  -18  have  multiple  intersections:  the  extra 
intersection*  are  labeled  G,  denoting  ghosts.  A crude  automatic  algorithm  has  been  used  in 
all  our  examples  to  distinguish  real  position  measurements  from  ghosts,  and  to  trace  the  com- 
puted path. 

2.  Azimuth  Krrors 

We  now  consider  the  effects  of  errors  in  azimuth  measurement.  Suppose  each  such  mea- 
surement <t>  is  replaced  by  a range  of  measurements  (i !>  - e,  <t>  t-  e) . where  e is  the  maximum 
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possible  error  in  4>.  Then,  each  possible  position  curve  is  replaced  by  a region  shaped  like 
a highway.  The  boundaries  of  this  region  are  just  the  possible  position  curve  rotated  by  *e 
and  - e degrees  around  the  sensor,  Figure  III  15-1  shows  these  possible  position  regions  for 
t 10,  IK,  and  6 with  z 0 and  e 1 ’ . Also  shown  are  the  computed  path,  plus  the  ghosts 
fort  18. 

[•'or  particular  t and  z.  the  intersections  of  the  S,  and  S , possible  position  regions  deter- 
mine an  intersection  region  where  the  plane  will  be  found.  This  region  generally  contains  all 
plane  locations  and  ghosts  found  in  the  previous  section  for  given  t and  z.  Basically,  ghosts 
are  the  result  of  nearly  tangentially  intersecting  possible  position  curves,  and  are  associated 
with  an  elongated  intersection  region  containing  the  true  position  plus  ghosts.  Noise  in  the 
azimuth  readings  can  cause  ghosts  anywhere  in  the  intersection  region. 

The  intersection  angle  of  the  and  possible  position  curves  indicates  whether  the  inter 
section  region  will  he  elongated,  for  mathematical  reasons,  we  consider  this  angle  to  have  a 
range  from  0 to  " , with  180“  being  the  usual  angle  for  tangential  intersection.  Angles  be- 
tween 180’  if  and  180“  i 2e  will  generally  result  in  elongated  intersection  regions,  whereas 
angles  near  00  or  i 70“  will  generally  result  in  small  rectangular  intersection  regions.  The 
intersection  angle  can  be  computed  as  follows.  From  the  point  at  which  the  plane  is  located, 
draw  three  velocity  vectors  (see  Fig.  Ill- 15-4): 

V Plane's  apparent  horizontal  velocity. 

Vj  Points  at  Sj  and  has  magnitude  (speed  of  sound)/c.os  , 

where  O,  is  the  elevation  angle  of  the  plane  at  the  given 
point  Idepends  on  z). 

v^  Similar  to  v( , but  for  and  0^. 

Then,  the  direction  of  the  Sj  possible  position  curve  going  through  the  point  is  Wj  v - vf . 
and  of  the  S2  curve  is  w^  "v  - v.,. 

As  an  aside,  we  indicate  why  the  above  procedure  works.  At  time  t,  the  plane  is  at  the 

given  position,  emitting  sound  received  by  Sj  at  time  t^  and  angle  0^.  At  time  t 4 dt,  the  plane 

has  moved  by  an  incremental  position  dt  57,  emitting  sound  received  by  at  time  t(  * dt^  and 

angle  01  + d 0(.  Bet  be  a vector  directed  along  the  possible  position  curve  for  Sf , such  that 

moving  an  incremental  position  dt  w(  from  the  original  position  will  also  correspond  to  the 

same  change  (dt^ , d ^ ) in  the  azimuth  curve  for  Sj . Consider  the  local  rectilinear  coordinate 

system  7 ,7,  ,7,  at  the  aircraft  location  associated  with  global  spherical  coordinates  with  origin 
J r <J>  »> 

S)(  where  Tr  is  the  unit  vector  directed  from  S,  to  the  given  position,  7^  is  the  unit  vector  in 
the  direction  of  increasing  <t>,  and  Tf)  is  the  unit  vector  in  the  direction  of  increasing  0.  Then, 
If  r^  is  the  horizontal  distance  between  Sj  and  the  plane  location, 

d0  j - dt(  v . T^)/r1 

<10,  - dtlWj  • T^l/rj 

dtj  = dt  + dt(  v • ir/speed  of  sound) 

dt,  = dt(  w,  ■ lr/speed  of  sound) 


These  last  equations  imply 


w.  ■ i . v ■ i. 
t <f>  <P 

w,  • i speed  of  sound  + v • i t 

I r ’ r 

The  equations  just  given  imply  Wj  v - v^,  where  Vj  is  chosen  so  that 
speed  of  sound 
0 

While  v is  the  apparent  horizontal  velocity  for  the  computed  path,  and  not  the  plane's  actual 
horizontal  velocity,  we  will  ignore  the  distinction  here  because  the  difference  between  these 
velocities  will  he  small  as  long  as  the  error  in  position  due  to  mis-estimated  height  z is  small, 
fir  the  position  were  accurate,  both  these  velocities  would  be  the  same,  because  they  give  the 
same  rates  of  azimuth  change  for  and  S^  l Note,  however,  that  an  actual  vertical  component 
to  the  plane  velocity  is  equivalent  to  a radial  horizontal  component,  because  it  projects  on  ir, 
but  that  this  equivalent  horizontal  component  is  different  for  each  sensor,  and  cannot  be  ex- 
pressed by  changing  the  apparent  horizontal  velocity  of  the  plane.  So,  our  calculations  are  not 
exact  unless  the  plane  actually  has  zero  vertical  velocity,  or  unless  the  assumed  height  z 0. 

For  given  z and  v,  curves  of  constant  intersection  angle  can  be  plotted  in  the  horizontal 
plane.  These  constant  intersection  angle  curves  have  sensors  as  end  points,  and  are  shown  in 
Fig.  III-H-I  for  z 0 and  the  actual  plane  velocity  V.  We  call  the  region  in  which  the  intersection 
angle  (ia)  lies  between  180  2e  and  180°  i 2e  the  region  of  linear  ambiguity,  because  in  that 
region  plane  position  can  only  be  localized  in  one  dimension. 

Comparison  of  Figs.  IJJ-H-3  and  I1I-H-5  suggests  that  the  intersection  angle  is  not  a very 
precise  measure  of  how  well  a plane  is  localized.  The  region  of  linear  ambiguity  is  nonetheless 
a useful  descriptive  concept. 

Figure  III-B-6  shows  the  computed  paths  of  the  four  comers  of  the  intersection  regions  for 
times  including  those  of  Fig.  Ill  -B-2.  The  intersection  regions  tend  to  lie  between  the  outermost 
of  these  four  paths.  One  clearly  sees  the  changes  in  length  of  the  intersection  regions  as  the 
intersection  angle  changes.  Near  the  region  of  linear  ambiguity,  the  corners  of  the  intersection 
regions  exchange  positions  relative  to  each  other,  with  their  computed  paths  crossing  over  each 
other. 

3.  Errors  in  Assumed  Height 

Figure  Ill-B-7  displays  the  computed  plane  paths  for  three  different  assumed  heights: 
z 0,  1 , and  2 km.  Not  much  change  is  observed  as  z changes,  except  in  the  region  of  linear 
ambiguity. 

Figure  II1-H-8  is  just  Fig.  III-B-2,  but  with  assumed  height  z ; 2 km  instead  of  z = 0.  Note 
that  not  only  have  the  ghosts  disappeared,  but  there  is  also  no  intersection  at  all  of  the  possible 
position  curves  for  t - -18  sec.  Errors  in  azimuth  can  also  eliminate  intersections. 

Possible  position  curves  for  t 0 and  z = 0,  t,  and  2 are  displayed  in  Fig.  I1I-B-9.  Changes 
in  z only  cause  changes  in  the  radial  distance  to  the  sensor,  and  therefore  change  the  possible 
position  curve  significantly  only  when  that  curve  is  perpendicular  to  the  radial  direction  from 
the  sensor.  Maximum  perpendicularity  is  achieved  by  fast  (nearly  Mach  1.0)  planes  flying 
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almost  toward  the  sensor.  Hut  if  such  a plane  is  between  the  two  sensors,  it  will  have  a nearly 
radial  possible  position  curve  for  the  other  sensor,  and,  consequently,  the  possible  position 
curve  for  the  other  sensor  will  not  change  much  with  changing  z.  So,  changing  z in  Fig.  Ill  - II -9 
moves  the  intersection  point  along  one  of  the  two  intersecting  possible  position  curves. 

( hanging  /.  induces  a radial  translation  in  a possible  position  curve  very  similar  to  the 
radial  translation  induced  by  changing  t.  Thus,  changing  z in  Fig.  Ill -11-9  moves  the  intersec- 
tion point  along  the  actual  path.  As  a result,  mis-estimating  z floes  not  mis-estimate  the  plane 
path,  hut  only  the  times  at  which  the  plane  is  at  various  points  on  that  path. 

For  planes  flying  between  the  sensors,  changing  z generally  results  in  a position  error 
parallel  to  the  plane  path. 

4.  Conclusions 

Time  series  azimuth  measurements  from  just.  two  sensors,  plus  an  estimate  of  plane 
height,  can  tic  used  to  locate  a plane  within  a box,  the  intersection  region  of  two  possible  posi- 
tion regains.  For  approximately  half  of  a typical  plane  path,  this  intersection  region  is  a nicely 
behaved  rectangular  box,  hut  elsewhere  the  region  is  elongated  and  may  become  U-shaped.  The 
major  exception  to  this  rule  occurs  when  the  plane  flies  along  the  straight  line  through  the  sen- 
sors, passing  very  near  both  sensors. 

For  two  sensors  10  km  apart,  and  planes  flying  no  higher  than  1 km,  the  effect  of  mis- 
estimating the  plane  height  is  usually  minimal. 
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Fig.  I1I-B-1.  Azimuth  angles  vs  time  measured  by  sensors  and  S, 
for  some  aircraft  to  be  located. 


Fig.III-B-2.  possible  position  curves  for  t from  - 10  to  —10  sec.  There 
is  one  curve  for  each  sensor.  Si  and  S^.  and  intersections  of  these  two 
curves  are  either  points  on  computed  plane  path  (nearly  straight  line)  or 
ghosts  (marked  G).  Assumed  plane  heightis  z 0,  actual  height  isO.l  km, 
plane  velocity  is  Mach  0.9. 


Fig.  IfI-B-3.  Possible  position  regions  and  intersection  regions  for  times  -30. 
-18,  and  -6  sec.  Computed  plane  path  (straight  line)  and  ghoBts  (marked  G)  for 
these  times  are  also  plotted.  Possible  position  regions  are  bounded  by  possible 
position  curves  with  assumed  errors  e of  +5°  and  -5°.  z 0 as  for  Fig.  I1I-B-2. 
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Fig.III-B-6.  Computed  paths  of  corners  of  intersection  regions  for  errors  e 
of  +5’  and  — 5 " in  azimuth  angle,  assumed  height  z - 0,  and  times  from —50 
to  450  sec  in  4-sec  increments.  Intersection  regions  generally  lie  between 
outermost  paths.  Kegion  of  linear  ambiguity  is  also  shown. 


Fig.  III-B-7.  Computed  plane  paths  for  assumed  heights  z = 0,  1,  and  2. 
Plane  location  is  determined  at  times  from  -50  to  +50  sec  in  4-sec  in- 
crements. Region  of  linear  ambiguity  is  also  displayed. 
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IV.  GENERA!.  SYSTEM  ANALYSIS 

One  fundamental  Issue  for  a surveillance  system  is  the  optimization  of  nodal  geometry. 
Should  a DSN  consist  of  a few  very  powerful  sensors  or  a large  numbi  r of  less-capable  sensors? 
During  this  reporting  period,  our  work  in  the  area  of  general  system  analysis  has  concentrated 
upon  a limited  aspect  of  this  question.  The  issue  is  a complex  one  involving  communications, 
system  cost,  performance,  reliability,  survivability,  and  other  factors.  An  initial  analysis 
has  been  completed  which  has  focused  upon  total  system  cost,  sensor  cost,  the  overlap  in  sen- 
sor coverage,  sensor  performance,  and  signal-path-loss  factors.  The  analysis  makes  explicit 
some  basic  system  trade-offs.  This  study  was  to  identify  and  clarify  issues  and  not  to  generate 
definitive  answers  to  specific  design  questions.  The  analysis  concentrated  upon  the  surveillance 
of  two-dimensional  regions. 

Suppose  an  area  A t 00  X 100  km  is  to  be  covered  using  sensors  which  each  have-  a target 
acquisition  range  of  R km.  Figure  IV-1  shows  the  number  of  sensor  sites  required  and  the  total 
system  cost  for  different  costs  per  sensor.  In  constructing  the  figure,  we  assumed  that  sensors 
can  cover  a disk  of  area  and  that  sensors  are  placed  on  a square  grid  of  points,  and  that  grid 
size  is  adjusted  to  the  largest  value  which  leaves  no  point  uncovered  by  at  least  one  sensor. 
Figure  IV-2  shows  this  arrangement.  Obviously,  other  packings  and  rules  of  coverage  could 
be  used,  and  we  have  selected  a reasonable  one  for  specific  discussion.  This  packing  gives  an 
area  coverage  redundancy  factor  p ir/2,  which  was  used  in  constructing  the  figure.  Commu- 
nication costs  are  assumed  included  in  the  sensor  site  costs.  Figure  IV-1  involves  all  basic 
system  attributes  of  interest  in  our  preliminary  analysis  of  sensor  performance  and  density 
issues  except  for  propagation  path  loss.  Some  discussion  of  this  figure  will  serve  to  motivate 
and  introduce  the  -subsequent  analysis. 

For  A and  p as  specified  above,  we  note  that  a system  using  sensor  sites  with  a range 
| of  2 3 km  and  costing  $500,000  each  would  result  in  a total  system  cost  of  5 million  dollars. 

Suppose  there  is  an  option  to  use  $50,000  sensor  sites.  If  the  range  of  the  less -expensive  sen- 
sor is  greater  than  7 km,  the  total  system  cost  will  be  less  for  the  large  number  of  less- 
expensive  sites.  This  issue,  trade-off  between  sensor  cost  and  range,  is  central  to  the  analysis 
which  follows.  The  analysis  is  to  clarify  issues  and  increase  understanding.  It  should  not  he 
interpreted  as  a definitive  solution  to  the  question  of  optimum  sensor  distribution. 
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A.  COS  T MODI; I . IN  TF1RMS  OK  SN  It  GAINS 

[•'or  a passive  sensor  system,  the  basic  signal-to-noiso  ratio  (SNIt)  at  a sensor  site  depends 
on  the  range  to  the  target  and  the  target  signal  strength.  The  sensor  output  SNIt  can  be  increased 
by  signal  processing,  reducing  system  noise,  the  use  of  multiple  sensors  combined  as  an  array, 
or  some  other  method  of  obtaining  directivity.  This  is  true  for  systems  involving  passive  sonar, 
atmospheric  acoustics,  lit  sensors,  seismic  sensors,  etc.  Kor  an  active  system,  the  basic 
SNIt  at  a site  depends  on  range,  the  power  radiated  in  the  direction  of  the  target,  and  the  target 
cross  section  which  defines  the  effectiveness  of  the  target  as  a re-radiator  of  power  in  the  di- 
rection of  the  receiver.  The  SNIt  can  be  increased  by  signal  processing,  increasing  the  power 
radiated  toward  the  target,  increasing  the  directivity  of  antennas,  or  reducing  receiver  noise 
levels.  Our  basic  sensor  cost  model  for  active  and  passive  systems  is  that  sensor  cost  is  a 
monotonic  increasing  function  of  the  SNIt  gain  of  the  sensor. 

Consider  a minimal  sensor  of  some  type  with  an  effective  operational  range  of  Rq  for  tar- 
gets of  interest.  The  value  of  lt(J  will  depend  upon  the  target  type,  but  we  assume  that  this  is 
fixed  for  our  analysis.  We  also  assume  that,  all  other  factors  being  fixed,  sensor  cost  is  a 
function  of  sensor  output  SNR.  For  this  purpose,  SNR  is  defined  as  a power  (squared  ampli- 
tudes or  mean  squared  values)  ratio.  Let  tl  be  the  SNR  gain  >f  an  enhanced  sensor  relative 
to  the  minimal  sensor.  The  model  for  the  cost  of  a sensor  is 


where  C.  is  the  incremental  cost  per  unit  of  SNR  gain,  and  c is  a positive  number.  Note  that 
the  reference  sensor  cost  is  Cq  + C^,  not  just  Cq.  In  what  follows,  we  will  deal  with  the  cases 
c t and  c = 1/Z. 

A few  examples  will  indicate  that  this  model  is  a reasonable  one  for  a number  of  situations 
i of  interest.  Consider  a passive  system  using  microphones,  seismometers,  or  hydrophones. 

I'he  most  common  way  to  improve  the  capability  of  an  individual  site  is  to  install  multiple  sen- 
sors and  to  o[)orate  the  site  as  an  array.  To  first  order,  it  is  reasonable  to  assume  that  the 
site  cost  will  be  a linear  function  of  the  number  of  instruments  used  in  the  array.  Rut  if  the 
signal  is  coherent  across  the  array  and  there  are  no  special  noise  conditions,  the  array  will 
effectively  reduce  the  background  noise  power  level  by  a factor  of  N — where  N is  the  number 
of  sensors  in  the  array.  A reduction  of  noise  power  by  a factor  of  N is  the  same  as  increasing 
the  output  SNR  by  a factor  of  G 1/N.  Thus,  our  cost  model  should  bo  applicable  to  a wide 
range  of  passive  systems  with  the  exponent  c set  equal  to  unity.  F'or  such  passive  systems, 
the  natural  reference  sensor  site  is  one  containing  a single  sensor  (microphone,  etc.). 

The  situation  is  a little  more  complicated  for  active  systems.  F’irst,  we  consider  situa- 
. tions  which  seem  reasonable  for  c - 1.  The  SNR  of  an  active  radar  or  sonar  can  be  increased 

by  increasing  the  peak  radiated  power  or  increasing  the  pulse  lengths  or  integration  times.  The 
SNR  is  linearly  proportional  to  peak  radiated  power  and  to  integration  time.  Increasing  integra- 
tion time,  particularly  in  a digital  system,  could  result  in  a linear  increase  in  processing  costs. 
This  would  tend  to  justify  our  cost  model  with  c unity.  It  is  also  reasonable  to  assume  that 
sensor  cost  will  increase  linearly  with  peak  power.  A simple  argument  to  support  this  is  to 
note  that  two  equal-power  sources  can  be  combined  to  give  an  effective  single  source  of  twice 
the  power  and  twice  the  cost  of  either  of  the  small  sources.  Thus,  an  increase  in  peak  |x>wer 
level  can  also  be  associated  with  a cost  model  having  c equal  to  unity.  Finally,  for  the  c = 1 
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case,  consider  that  a cylindrical  array  of  horns  is  used  for  either  transmitting  or  receiving, 
hut  not  both.  To  first  order,  the  antenna  power  gain  will  be  directly  proportional  to  the  number 
of  horns  needed  to  cover  all  azimuths,  and  so  will  the  cost  of  such  a configuration.  If  horns  are 
used  for  both  transmission  and  reception,  the  cost  still  remains  a linear  function  of  the  number 
of  horns  but  the  SNR  gain  goes  as  the  square  of  that  number.  This  is  the  important  case  for 
which  c = t/2. 

Our  model  of  sensor  site  cost  assumes  that  communication  costs  are  included  in  (,‘0.  This 
is  reasonable  since  we  assume  communication  by  a distributed  packet  switched  communication 
system  of  the  Packet  Radio  Net  (PRN)  class  and  tnat  the  !*RN  single-hop  range  is  at  least  as 
large  as  the  distances  between  adjacent  sensor  sites.  For  such  a system,  C(J  is  bounded  from 
below  by  the  cost  of  a PRN  site.  Systems  using  hardwire  communication  or  needing  extra  re- 
peaters would  require  more-complicated  modeling,  and  would  generally  increase  the  pressures 
to  select  more-powerful  sensor  sites. 

With  the  above  model  for  single  sensor  site  cost,  the  model  for  total  system  cost  is  very 
simple.  Specifically,  as  used  in  our  introductory  figure, 

CT  = pCA/irR2 

where  A is  the  area  to  be  covered,  irR2  is  the  area  of  the  disk  covered  by  a single  site  with 
effective  range  R,  and  p is  a redundancy  factor. 
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H.  CONVERSION  OF  COSTS  TO  A RANGE  FUNCTION 


i 


If  a power  gain  of  G > t Ih  obtained  for  a sensor  site,  this  will  increase  the  effective  range 
from  ItQ  to  some  K > ly  The  amount  of  increase  in  range  depends  on  the  path  loss  experienced 
by  the  signals.  A convenient  model  for  path  loss  is  that  the  power  loss  as  a function  of  range  is 
of  the  form 

UK)  K* 

where  A is  a positive  real  number.  For  this  model,  the  extended  range  of  a sensor  as  a result 
of  a gain  G in  SNK  is  the  solution  of 

G = m/R0ix  . 

Using  this,  the  cost  of  a sensor  site  becomes 

c c^  + eyit/iy" 

where  n Ac  depends  upon  both  the  path  lose  model  exponent  A and  the  SNK  gain  cost  exponent  c. 
With  this  the  total  system  cost,  in  terms  of  the  range  of  a single  tensor,  becomes 

Ct(K)  (pA/ir)  |C0R'2  + C1Kn‘2K0"n|  . 

In  the  next  section,  we  discuss  the  mathematically  optimum  selection  of  a sensor  range. 

To  facilitate  that  discussion,  we  include  here  a brief  review  of  the  physical  significance  of  our 
path  loss  model  for  different  exponents.  Imagine  a target  radiating  energy  from  a point  in  a 
homogeneous  three-dimensional  space.  Such  a situation  involves  spherical  spreading.  Signal 
power  density  is  an  inverse  squared  function  of  the  distance  from  the  radiator.  This  follows 
directly  from  the  fact  that  the  surface  area  of  a sphere  is  4»K^.  Thus,  our  path  loss  model 
with  A 2 is  appropriate  for  a passive  system  with  the  signal  subject  to  s|>herical  spreading. 
Other  passive  sensor  situations  lead  to  less  path  loss.  For  example,  acoustic  signals  in  the 
atmosphere  or  ocean  can  be  trapped  within  a two-dimensional  waveguide  so  that  the  spreading 
of  signals  is  in  two  dimensions  rather  than  three.  This  is  the  case  of  cylindrical  spreading,  anil 
corresfionds  to  a path  loss  model  with  A = i.  Finally,  one  can  have  more  path  loss  than  spher- 
ical spreading  so  that  A > 2.  This  can  happen,  for  example,  if  the  medium  itself  is  very  lossy, 
if  variations  in  the  medium  result  in  strong  shadow  zones,  or  if  there  are  other  physical  ob- 
structions to  the  propagation.  A nominal  value  of  A for  scattered  or  diffracted  signals,  which 
are  of  some  interest  to  us,  is  A = 4. 

The  situation  with  active  systems  is  similar,  but  we  must  use  path  loss  models  with  twice 
the  exponent  value  used  for  passive  systems.  For  example,  consider  a radar  system.  The 
physical  propagation  involved  is  spherical,  exclusive  of  terrain,  foliage  loss,  etc.  However, 
the  spherical  spreading  occurs  for  the  signal  radiated  from  the  transmitter  and  it  also  holds 
for  the  signal  reradiated  from  the  target.  Thus,  radar  Involves  squared  spherical  spreading 
and  corresponds  to  the  path  loss  model  A = 4.  In  the  case  that  the  radar  is  line-of-sight  (LOS) 
limited  beyond  ly  an  appropriate  value  for  A would  be  «o.  The  value  for  diffracted  radar  would 
be  A = 8. 
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c. 


OPTIMUM  STKAIKGY  AM)  DISCISSION 


Given  lt0.  C0,  Cj,  and  c.  we  want  to  determine  if  a particular  option  for  enhancing  a 
sensor  should  be  used.  Clearly,  if  C.(  (III  is  less  than  ('  ( (lt()l  for  some  It  lt(),  then  the  site 
should  be  upgraded.  Now  we  note  that  if  n Ac  is  less  than  or  equal  to  i,  then  C.p(K)  is  a 

monotonie  decreasing  function  of  |(.  In  this  situation  tie  site  .hould  definitely  be  enhanced. 

Of  n :rsc.  this  strategy  is  valid  only  as  long  as  the  model  is  valid.  If  n is  greater  than  2, 
then  the  value  of  It  which  minimizes  the  total  sy  tern  cost  is  obtained  by  setting  the  derivative 
of  the  cost  to  zero.  This  gives 

It*  K0  l^'oA't'n  ~2,|l/n 

Depending  u|mn  the  value  of  £('  /('  ( n 2).  tin  strategy  may  be  to  enhance  the  site  or  reduce 
its  eapahilitie  . l arge  valines  of  n will  tend  to  push  toward  small  sites,  while  small  values 
of  n will  have  the  opposite  effect.  In  terms  of  propagation  loss,  a similar  trend  is  seen.  That 
is.  large  propagation  loss  mitigates  toward  many  small  sensor  sites,  and  small  propagation 
loss  favors  a few  powerful  sites.  Note  that  the  redundancy  factor  p does  not  enter  into  the  op- 
timum choice  of  K,  although  it  clearly  affects  system  costs  and  performance. 

Table  IV-1  gives  a more-physical  interpretation  of  the  implications  of  minimizing  our  ex- 
pression for  system  cost  for  a number  of  different  circumstances.  In  the  table,  the  use  of 
directivity  for  both  transmission  and  reception  is  the  equivalent  of  using  an  SNK  cost  model 
exponent  c - 1/2.  All  other  SNK  improvement  strategies  correspond  to  c 1. 

First  note  that  the  use  of  directional  antennas  for  both  transmission  and  reception,  within 
the  context  of  our  model  and  analysis,  is  a very  cost-effective  thing  to  do.  For  this  SNK  im- 
provement option,  one  continues  to  increase  directionality  as  much  as  possible  so  long  as 
cost  is  linear  with  directionality  and  transmission  loss  is  no  worse  than  spherical.  For  radar, 
the  obvious  factor  which  impacts  both  of  these  is  the  LOS  limit.  When  this  limit  is  reached,  it 
can  be  interpreted  as  a sudden  switch  to  a propagation  loss  model  which  is  much  worse  than 
spherical.  This  will  cause  us  to  stop  trying  to  increase  range  by  directionality.  Also,  one  way 
to  improve  the  LOS  situation  is  to  use  well-sited  and  elevated  radars.  This,  in  turn,  may  cause 
the  radar  cost  to  increase  at  a much  greater  rate  than  that  which  is  in  our  model.  In  any  case, 
for  active  systems  which  can  have  directional  gain  increased  at  linear  cost,  we  see  that  the  trend 
is  toward  larger  sensors  up  until  other  secondary  factors  make  further  enhancement  practical  or 
desirable. 

Now  consider  the  situations  in  which  the  SNK  improvement  strategy  results  in  a linear  re- 
lationship between  SNK  gain  and  sensor  cost.  That  is,  consider  all  SNK  improvement  strategies 
except  using  directivity  for  both  transmission  and  reception.  This  includes  all  passive  systems 
and  some  of  the  active  options.  Also  exclude  the  case  of  LOS  limited  propagation  since  it  Is 
clearly  a special  case  insensitive  to  SNK  improvement  strategy. 

For  active  sonar  the  boundary  where  we  always  choose  to  enhance  the  sensor  moves  from 
the  spherical  to  cylindrical  propagation  loss  regime.  That  is,  if  loss  is  greater  than  cylindrical, 
then  a limit  to  sensor  enhancement  will  be  reached.  We  point  this  out  since  cylindrical  propaga- 
tion is  an  important  mode  for  underwater  sound  which  is  often  trapped  in  a channel. 

For  passive  systems  in  general,  the  boundary  situation  is  the  case  of  spherical  spreading. 
This  is  true  for  passive  sonar  systems  and  for  atmospheric  acoustic  systems.  For  example, 
if  sound  is  trapped  in  a channel  resulting  in  cylindrical  spreading,  then  Individual  sensor  sites 
clearly  should  be  enhanced. 

i 
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TABLE  IV-1 


OPTIMAL  SENSOR  RANGES  FOR  VARIOUS  SYSTEMS,  SNR  IMPROVEMENT 
STRATEGIES,  AND  PROPAGATION  LOSS  MODELS 


System 

SNR 

1 Improvement 
Strategy 7 

Propagation 

Loss  Model 

Cost 

Exponent, 

n 

- 

Optimal  Sensor 
Range,  R 

RADAR 

D 

Spherical  or  less 

2 or  less 

Increase  R 

RADAR 

D 

More  than  spherical 

>2 

R = R0l2C0/C1(n-2))’/n 

RADAR 

Any 

LOS  limited 

JC 

R = R0 

Active  SONAR 

D 

Spherical  or  less 

2 or  less 

Increase  R 

Active  SONAR 

D 

More  than  spherical 

>2 

R=  R0l2C0/C](n-2)|'/n 

Active  SONAR 

A,  B,  C 

Cylindrical  or  less 

2 or  less 

Increase  R 

RADAR 

A,  B,  C 

Spherical 

4 

R = R0(C0/C,)'/4 

Active  SONAR 

A,  B,  C 

Spherical 

4 

R=R0(C0/C,),/4 

Passive 

A,  C 

Spherical  or  less 

2 or  less 

Increase  R 

Passive 

A,  C 

More  than  spherical 

>2 

R = R0l2C0/C1(n-2)i1/n 

Passive 

A,  C 

Scattered 

4 

R " Ro(Co/ci)1/4 

Passive 

Any 

LOS  limited 

JO 

l 

o 

O' 

II 

Cd. 

t SNR  Improvement  Strategy  Codes: 

A.  Directivity  on  transmission  or  reception  only.  In  the  case  of  arrays, 
increase  the  number  of  array  elements. 

B.  Increase  transmitter  power  of  active  system. 

C.  Increase  pulse  lengths  or  integration  time. 

D.  Directivity  used  on  both  transmission  and  reception. 


Finally,  Table  J V - 1 includes  three  situations  wjth  n 4.  The  passive  system  case  with 
scattering  ns  the  principal  propagation  in  of  particular  intercut  to  us  since  it  may  be  a reason- 
able model  for  acoustic  surveillance  of  low-flying  aircraft  once  tin-  geometric  shadow  zone  is 
entered.  In  this  case,  we  see  that  our  model  asserts  that  sensor  performance  should  be  en- 
hanced as  long  as  ('q  is  greater  than  C,.  The  optimum  sensor  is  that  for  which  they  are  equal. 
One  interpretation  of  this  is  as  follows.  In  tin-  model,  C.q  includes  the  irreducible  cost  of  a 
sensor  site.  For  example,  it  includes  the  cost  of  the  required  communication  system.  The  sen- 
sor itself  may  add  to  C.q,  hut  it  is  clear  that  an  optimum  sensor  site  should  cost  at  least  twice 
as  much  as  that  irreducible  cost. 
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